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A  progress  report  "Insulated  Test  Road  ...  State  Road  26"  by  Richard  P. 
Stulgis,  Graduate  Assistant  in  Research  on  our  staff,  is  presented  to  the 
Board  for  acceptance.  This  research  was  directed  by  Professor  C.  W,  Lovell, 
Jr.  of  our  staff. 

A  plan  has  been  developed  for  the  construction  of  a  thermally  insulated 
test  road  situated  on  a  relocation  of  Stete  Road  26,  Just  west  of  the  town 
of  Rossville.  The  proposed  installation  has  three  sections  (two  insulated 
and  one  conventional  design  for  control);  each  section  being  200  feet  in 
length.  The  foamed  plastic  thermal  insulation  is  to  be  used  in  a  1-inch 
tbicimess  at  the  base  of  a  normal  section,  and  in  a  lg~inch  thickness  in  a 
reduced-thickness  section.  The  relative  effectiveness  of  the  insulation 
in  preventing  frost  penetration  through  the  pavement  sections  will  be  evalu- 
ated by  monitoring  105  subsurface  thermistors.  Surface  deflection  measure- 
ments are  also  recommended. 

This  report  is  in  fulfillment  of  a  study  proposal  approved  by  the 
Board  on  September  27,  1967^  The  test  road  is  planned  for  construction 
before  the  Winter  of  1968-69,  Detailed  plans  for  collecting  data  and 
evaluating  performance  axe   to  be  worked  out  Jointly  by  selected  personnel 
of  the  Research  and  Training  Center  and  the  Joint  Highway  Research  Project, 
and  several  additional  reports  to  the  Board  are  anticipated  over  the  next 
several  years. 
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ASSTBACT 

'Qx's  requisite  planning  for  an  experimental  section  of  Indiana 
hi^way,  ^Incorporating  fosun-plastic  Insulating  layers  to  attenuate 
fiost  penetraticmj  •was  developed  o 

Vftth  the  help  of  Indiana  state  hi^vray  officials,  a  flexible 
pavement  construction  project  (Project  Noo  S2Yl)  on  State  Road  26 
was  selected^  vithin  MHaldh  the  proposed  field  installation  vould 
he  located o  Two  insulated  sections  and  one  non°insuXated  section 
coogprisa  the  test  iastallationo 

Selection  of  the  insulation  thiclsnesses  to  be  used  "was  based 
oa  the  results  of  a  coa^xxter  analysis  developed  at  Purdue  IfaiversitjTo 
Eiis  program^,  under  the  condition  of  one-dimensiocal  heat  flow  by 
conduction^  pred3.cts  the  distribution  of  tengperature  ■with  time 
throu^out  a  layered  losdluittc.  Input  data  utilized  were  based  on 
climatic  records  for  the  area  in  ^ich  the  site  is  located  and  on 
estimated  thernal  propex'ties  of  the  cras^jonent  layers  of  Idae  proposed 
hi^'ijay  cross-'Sectiono 

Methods  of  evaluating  the  thermal  and  structural  performance 
of  the  insulated  and  non-insulated  sections  irere  recommended o  Hxermal 
perfo2iEance  is  to  be  evaluated  by  means  of  thermistors  strategically 
placed  in  each  of  the  three  sections  and  structural  performance  by 
aieans  of  Benkelman  beam  tests  <> 

Special  construction  procedures  "were  recomraended  for  the  field 
instaUAtion^  due  to  the  ejcperiioental  nature  of  the  project  o 


msRowcTzm 

de  use  of  a  thermal  barrier  benea-th  a  transportation  route; 
in  oixter  to  minimize  ftrost  damagej  has  received  a  recent  linpettjs 
with  the  availability  of  foamed  plastics  vihich  are  excellent 
insulators  o 

Research  and  development  effort  over  the  past  decade  has 
ahovm  that  the  use  of  such  materials  is  technolo^cal  ly  feasible  and 
that  such  use  can  be  economically  desirable  o  As  market  applications 
Increase,  iBiproveisents  in  the  foamed  plastic  properties  and  in  the 
technology  of  its  placement  can  be  e3q>ected  ooo  sis  veil  as  some 
reduction  in  its  as-placed  costo 

A  number  of  hifi^n^y  departments  are  in  the  process  of  gaining 
first  hand  experience  '.jith  the  insulated  design  concept  by  means 
of  test  roads  o  IMs  is  probably  a  necessary  step  for  each 
department,  since  soils,  loads,  environment  and  service  reqvilremeats 
vaiy  considerably  with  and  within  the  various  political  units  o 


PURPOSE 

The  purpose  of  the  study  was  to  develop  the  requisite  planning  for 
an  e:q)erimeatal  section  of  Indiana  highway  incorporating  foam-plastic 
insulating  layers,  including  appropriate  testing  and  evaluating  procedures. 

Although  field  test  installations  have  already  been  made  in  the 
northern  Uaited  States,  Cauada,  and  in  Europe, and  have  demonstrated  good 
performance,  the  Indiana  installation  has  two  definite  objectives.  The 
first  is  to  obtain  first-hand  experience  with  insulated  constiruction  and 
performance  of  the  insulated  design.  Secondly,  the  rather  complete  in- 
strumentation of  the  project  will  permit  both  validation  and  extension  of 
extant  analytical  models  for  thermal  pavement  design. 
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LOCATION  AND  DESCRIPTION  OF 
FIELD  TEST  INSTALLATION 

The  proposed  experimental  test  section  of  roadway  will  be  constructed 
in  a  cut  area  located  within  the  3.1  mile  length  of  flexible  pavement  con- 
struction proposed  on  State  Road  26  (S-Itoject  No.  271).  As  shov/n  on  the 
project  locus  (Figure  l),  this  construction  job  is  located  in  northwestern 
Clinton  County  near  Rossville,  Indiana. 

The  project  lies  in  the  north  central  portion  of  the  T5.pton  Till  Plain 
Physiographic  Subsection  of  Indiana  (^).  Rolling  and  hilly  topographic 
features  found  in  the  project  area  are  attributed  to  the  proximity  of  a 
se^ient  of  the  Bloomington  Morainic  system  which  is  located  to  the  south 
of  State  Road  26  in  an  east  and  west  direction.  An  exploratory  program, 
consisting  of  hand -borings,  carried  out  by  the  Division  of  Materials  & 
Tests  of  the  Indiana  State  Highway  Commission  (ll)  found  the  subgrade 
soils  in  the  vicinity  of  the  test  installation  to  be  of  the  A»6  group 
(AASHO  Classification).  The  textural  classification  of  these  soils  ranges 
from  silty  clay  to  clay  and  thus  a  frost  susceptible  soil  appears  to  com- 
prise the  subgrade,  at  least  to  the  depths  to  which  the  hand  borings  were 
carried . 

This  location  was  selected  on  the  basis  of  the  following  factors: 

(1)  Anticipated  subgrade  soils  in  this  area  are  classified 
as  Fk   (Corps  of  Engineers  Classification  System),  whidi 
is  highly  frost  susceptible. 

(2)  The  most  severe  test  for  the  insulation  would  be  in  a  cut 
area  where  the  soil  would  have  easier  access  to  water. 
The  proposed  construction  in  this  location  would  entail 
cuts  of  from  one  to  five  feet. 


1.  Numerals  underlined  in  parentheses  refer  to  entries  in  the 
Bibliography,  jage  51. 


(3)  The  road  is  currently  posted  as  "Rovigh  Pavement"  and  it 

has  been  determined  that  this  deterioration  is,  in  part, 

caused  by  ftost  actiono 
(U)  The  location  is  close  to  Lafayette,  which  will  facilitate 

trips  to  the  site  to  retrieve  data, 
(5)  Construction  of  the  test  installation  will  begin  in  the 

Summer  of  I968. 
As  shown  on  the  plan  (Figure  2),  the  experimental  roadway  will  be 
constructed  between  Station  100  +  00  and  Station  I06  +  00  of  the  con- 
tract section^  The  test  axea  will  be  divided  into  two  teat  sections  and 
one  control  section,  each  200  feet  in  lengrfch.  The  test  sections  are  de>" 
signated  "A"  and  "B"c  Section  "A"  will  include  a  one  inch  thick  insulating 
layer  and.  Section  "B"  will  include  a  one  and   one  half  inch  thick  insulating 
layer,  the  insulating  layer  in  both  cases  being  placed  on  the  subgradeo  The 
depth  of  this  layer  below  finished  grade  \f;ill  be  20  inches  in  Section  '"A" 
(which  will  utilize  the  normal  pavement  section  for  the  project)  and  Ik 
inches  in  Section  "b"  (which  will  elimin0.te  the  six  inch  subbase  layer)  o 
The  width  of  the  insulating  layer  will  be  3!*  feet  in  Section  "A"  and  1^6 
feet  in  Section  "B".  The  control  section, utilizing  the  nDrma;L  design 
for  the  project, will  be  located  between  the  ttro  test  sectioa';^  and  is  desig- 
nated Section  "C"o  Transitions  between  the  insulated  sections  and  the  con- 
trol and  normal  roadway  sections  will  have  to  be  provided,  diie  to  differ- 
ences in  subgrade  elevations  between  these  sections »  Recommended  transitions 
are  as  shown  in  Figure  3°  Vertical  strings  of  thermistors  will  be  installed 
throughout  the  cross  section  at  the  center  of  each  of  the  tliree  sections, 
transverse  to  the  roadway  center  lineo 
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HESIGN  PROCEDUEE 

The  design  of  the  field  test  installation  consisted  of  determining 
the  number,  type  and  length  of  pavement  sections,  the  thickness  and  depth 
of  the  insulating  layer,  the  width  of  the  insulation,  the  location  of  tem- 
perature sensing  elements,  and  other  subordiEate  details. 

numbers  Type  and  Length  of  Sections 

It  is  recommended  that  the  installation  be  coarprised  of  three  parts, 
two  insulated  sections  (one  utilizing  the  normal  pavement  section  and  the 
other  elimisaaticg  the  six  inch  subbase  layer)  and  one  control  or  normal 
design  section. 

Section  "A"  affords  a  direct  comparison  between  the  perfoimance  of 
equal  thicknesses  of  an  insulated  section  and  a  non-insulated  section 
(Control  Section).  Section  "B"  is  reduced  in  thickness  relative  to  the 
normal  section  by  omission  of  the  six  inch  subbase. 

The  insulation  serves  two  purposes?  (l)  It  prevents  fSost  penetration 
into  frost  susceptible  subgxade  soils,  thereby  preventing  frost  heaving  and 
weakening  of  the  pavonent  upon  thawing;  and  (2)  by  avoiding  a  reduced 
strength  condition  in  the  sub©:ade,  a  thinner  pavement  may  be  usedo  The 
cost  of  the  insulation  (placed)  must  be  offset  by  savings  in  aggregate 
and/or  im^oved  pavement  performance,  in  order  to  make  the  insulation 
method  economically  competitive  with  conventional  methods ,  Section  "B" 
is  infcended  to  demonstrate  the  reduction-in- thickness  concept,  although  it 
is  unlikely  that  said  reduction  is  sufficient  to  produce  a  net  saving  in 
construction  cost.  In  reducing  the  pavement  thickness,  however,  twc  ques- 
tions arose.  The  first  concerned  the  structural  adequacy  of  the  reduced 


thickness  pavement  section,  and  the  second,  the  stress  level  transferred 
to  the  insulating  layer  due  to  wheel  loads.,  A  CBR  method  of  analysis  (39),  for 
the  anticipated  soil  and  loading  conditions,  attested  to  the  structural 
adequacy  of  the  reduced  section  a  With  r«?'5p«>ct  to  the  second  question, 
advice  was  sought  from  a  representative  of  the  major  manufacturer  of 
highway  insulation  i^.     Previous  experience  indicates  a  minimum  flexible 
pavement  section  consisting  of  three  inches  of  bituminous  material  and 
twelve  inches  of  granular  base  will  hold  the  vertical  stress  on  the  insu- 
lation to  a  desired  limit  of  15  pslo  The  proposed  Section  "B"  is  three 
inches  deficient  with  respect  to  the  twelve  inch  gravel  minimum  but  exceeds 
the  lower  Irlmit  of  three  inches  of  bituminous  material  by  two  inches «  It  was 
concluded  tliat  the  reduced  section  would  provide  the  minimum  protection  de- 
sirable with  respect  to  the  stress  level  transferred  to  the  insulating  layerc 

Teat  sections  must  be  long  enough  to  relegate  end  effects  to  a  minor 
role„  Figure  k   compares  the  thermal  conditions  which  would  exist  in  a  con- 
tinuous section  with  those  of  a  relatively  short  test  section c 

Based  on  previous  research  C^),  a  test  section  50  feet  long  is  adequate 
for  thermal  evaluations.  However,  longer  sections  are  required  for  pavement 
performance  ratings »  Therefore,  within  this  framework  and  the  practical 
constraints  of  the  project,  a  length  of  200  feet  was  chosen  for  each  section » 
Also,  a  total  project  length  of  60O  feet  allows  the  entire  installation  to 
be  contained  within  a  single  continuous  cut,  which  is  slightly  greater  than 
600  feet  in  length.  Thus,  non-uniformity  of  subgrade  soil  conditions,  is 
held  to  a  practical  minimum , 

It  should  also  be  pointed  out  that  normal  construction  procedure  in 
this  cut  could  Involve  a  subgrade  treatment  of  undercutting,  removal  and 
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PLAN     VIEW 


A  CONTINUOUS    HIGHWAY 
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replacement  of  certain  fi-ost  susceptible  subgrade  material,  indicating  an 
additional  economic  factor  in  favor  of  the  insulation. 

Thickness  and  Depth  of  Insulating  Layer 
The  insulating  layer  should  be  placed  directly  on  the  subgrade o  If 
granular  materials  are  placed  between  the  insulating  layer  and  the  subgrade, 
heat  losses  in  the  lateral  direction  increase.  Conversely,  placing  the 
insulating  layer  on  a  subgrade  which  has  a  relatively  high  heat  capacity, 
reduces  heat  losses  in  the  lateral  direction  and  one~dimensional  heat  flow 
is  apt  to  constitute  a  reasonable  model  of  the  actual  phenomenon.  Placing 
the  insulating  layer  "deep"  in  the  pavement  section  is  desirable  not  only 
from  the  thermal  standpoint,  but  also  structurally  since  stresses  on  this 
layer  are  minimized  {h)» 

To  estimate  the  thickness  of  insulation  required  in  the  test  sections, 
a  computerized  solution,  developed  principally  at  Purdue  University,  was 
used  C2)»  This  program,  under  the  condition  of  one-dimensional  heat  flow 
by  conduction,  predicts  the  distribution  of  temperatures  throughout  a 
layered  medium  with  time,  by  means  of  a  finite  difference  technique »  There 
are  essentially  no  limitatioos  on  the  functional  form  of  initial  and  boundary 
conditions,  or  on  the  variation  of  physical  and  thermal  properties  of  the 
layered  system »  Pertinent  parts  of  the  program  and  an  exaniple  of  its  appli- 
cation to  Section  "A"  are  included  in  Appendix  Ao 

The  predictive  capability  of  the  program  is  taapered  by  the  quality  of 
the  input,  which  includes  the  dimensions  of  the  cross  section,  the  physical 
(and  thermal)  properties  of  the  layers,  the  air-surface  transfer  factor, 
an  air  temperature- time  function,  temperature  versus  time  at  a  relatively 
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great  depth,  and  an  initial  temperature -depth  function.  See  Figure  27 
of  Appendix  A.  The  procedures  used  in  selecting  the  various  input  data 
and,  where  applicable,  an  assessment  of  the  effect  of  the  chosen  values 
on  the  predicted  temperature  distribution  with  time,  are  described  beloi*, 
with  specific  reference  to  Section  "A". 

Dimensions  of  the  Cross  Section.  The  starting  point  is  the  normal  design 
section.  A  value  Judgement  is  made  relative  to  position  and  thickness 
of  the  insulation  layer  and  a  reduction  in  thickness  of  granular  material 
in  the  normal  section.  Having  done  this,  the  components  of  the  pavement 
section  are  further  subdivided  into  any  number  of  incremental  thicknesses 
desired.  The  variation  of  temperature  with  time  at  the  center  of  these 
increments  is  computed.  For  example,  in  Figure  27,  selection  of  increments 
five  and  seven  will  yield  the  temperature  one  inch  above  and  one  inch  below 
the  insulation.  By  repeating  the  sequence  of  data  cards,  several  problems 
can  be  run  in  one  batch  and  thus,  several  thicknesses  of  insulation  may 
be  investigated  at  one  time,  A  thermally  acceptible  section  is  one  which 
essentially  prevents  penetration  of  the  32  F  isotherm  through  the  insulation. 
Structiaral  and  econoniic  factors  must  also  be  judged. 

Physical  (and  Thermal)  Properties  for  Each  Layer.  Values  for  the  unit 
weight,  initial  water  content,  volumetric  heat  and  thermal  conductivity 
were  selected  for  each  of  the  component  layers  of  the  pavement  section. 
Values  for  the  bitmainous  surface  and  base,  compacted  aggregate  base, 
and  subbase  were  selected  from  values  for  similar  materials  used  in  the 
state  of  Maine's  test  installation.  Values  for  the  insulation  were  based 
on  reported  values  for  Styrofoam  HI  Plastic  Foam,  manufactured  by  The 
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Dow  Chemical  Company  (8),  v/hile  values  for  the  subgrade  were  based  on 
reported  values  for  Crosby  silty  clay  (17).  With  the  exception  of  the 
insulation,  the  selected  values  are  rather  approximate.  The  predictions 
may  be  examined  as  to  sensitivity  to  vaxiation  in  any  factor.  For  example, 
Figure  5  illustrates  the  effect  of  varying  water  content  in  the  base. 

Air-Surface  Transfer  Factor.  This  factor  attempts  to  relate  pavement 
surface  temperature  to  air  temperature.  This  relation  is  coo^lex,  being 
influenced  by  transient  local  environmental  conditions.  It  is  expressed 

as  a  ratio  of  (  ^airSmStSe"^"^ '    ^^®  ^^*^°  ^^^®  assumed  for  this 
analysis  is  0.99 .  ^^   to  the  fact  that  the  pavement  surface  will  gener- 
ally exist  at  a  temperature  which  is  higher  than  the  mean  air  temperature, 
the  selection  of  0.99  ^ox   the  air-surface  transfer  factor  is  consearvative . 

Initial  Conditions.  In  actuality,  the  computer  solution  generates  a  three- 
dimensional  surface,  on  which  any  point  represents  a  temperature  at  a 
given  depth  and  time .  From  Figure  27,  it  can  be  sia-mised  that  although 
the  upper  and  loiter  boundary  conditions  V7ill  greatly  affect  the  shape  of 
the  surface  generated,  the  initial  conditions  will  affect  the  functional 
relationship  only  at  short  times,  i.e.,  with  time  the  effects  of  the 
initial  conditions  greatly  diminish.  Therefore,  although  it  is  desirable 
to  predict  the  initial  temperature  distribution  throughout  the  pavement 
section  as  close  to  the  actual  as  practicable,  the  error  associated  vrith 
any  reasonable  assumption  should  not  be  significant.  Thus,  it  ^vras  assumed 
that  the  initial  temperature  distribution  was  constant  with  depth,  at  a 
value  of  50°F. 


1.  Including  a  base  coxirse  and  a  subbase  course. 
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Upper  Boundary  Coaditions.  In  the  absence  of  air  teinperature  data  for  the 
actual  test  site,  the  upper  boundary  condtions  were  determined  from  mean 
daily  air  temperature  data  taken  at  West  Lafayette  (Purdue  University). 
It  was  decided  to  use  the  freezing  index  as  a  measure  of  the  severity  of 
specific  X7inters  to  aid  in  choosing,  for  design,  the  most  severe  winter. 
After  collecting  temperature  data  and  computing  values  of  freezing  index 
for  the  winters  of  the  past  ten  years,  the  Winter  of  I962-I963  was  selected 
as  the  most  critical  one,  having  a  freezing  index  of  127^.0  degree-days. 
Figure  6  shows  the  curve  of  cumulative  degree-days  for  the  Winter  of  I962- 
1963,  and  Appendix  B  contains  the  computer  program  used  in  computing  the 
values  of  freezing  index.  Also  included  in  Appendix  B  is  the  Output  for 
the  Winter  of  1962-I963,  from  which  the  curve  in  Figure  6  was  constructed. 
The  freezing  index  is  an  admittedly  approxixoate  way  of  selecting  the 
critical  upper  boundary  condition,  although  it  is  one  that  can  be  rather 
simply  applied  for  any  location  of  air  temperature  record.  It  is  possible 
to  have  a  deeper  frost  penetration  with  a  slightly  lower  freezing  index, 
if  the  duration  of  the  freezing  period  is  longer.  However,  in  this  analy- 
sis, the  freezing  index  of  the  Winter  of  I962-I963  was  i»00  degree-days 
greater  than  the  v/inter  x^ith  the  next  highest  freezing  index,  the  duration 
of  the  freezing  period  being  about  the  same  for  both  years.  Accordingly, 
the  selected  value  is  probably  a  conservative  one. 

Iiower  Boundary  Conditions.   The  influence  of  the  lo^rer  boundary  is  illustrated 

by  example  in  Figure  7.  This  condition  is  a  difficult  one  accurately  to  assess, 

and  the  instrumentation  recommended  in  this  report  should  provide  valuable 

1.  The  freezing  index  is  defined  as  the  difference  between  the  maximum  and 
minimum  point  on  the  curve  of  cumulative  degree-days,  a  degree-day  represent- 
ing one  day  with  a  mean  air  temperature  one  degree  above  or  below  freezing. 
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insight  with  respect  to  it.  For  the  purposes  of  this  analysis,  the  lower 
boundary  was  taken  as  a  constant  tenrperature  value  of  yTT  at  a  depth  of 
eight  feet  below  the  roadway  surface.   A  value  of  56°F  at  a  depth  of  100 
inches  below  the  pavement  surface  was  used  for  the  lower  boundary  in  pre- 
dicting the  performance  of  Maine's  test  installation  (15) >  and  correlation 
between  the  predicted,  and  the  actual  values  of  temperature  was  quite  good« 
Therefore,  Indiana  being  located  further  south  than  Maine,  it  would  appear 
that  a  value  of  50  F  is  suitable » 

Before  presenting  the  results  of  the  analysis,  it  should  ba  pointed 
out  that  despite  the  error  that  can  be  introduced  through  the  assuniptions 
made  in  selecting  input  data,  good  correlation  between  predicted  and  actual 
temperature  distributions,  such  as  shown  in  Figures  8  and  9,  can  be  obtained o 

The  results  of  the  analysis  performed  for  Sections  "A"  and  "B"  indicate 
that  a  thickness  of  one  inch  and  one  and  one-half  inches,  respectively,  will 
provide  adequate  thennal  protiection.  Figures  10  and  11  depict  the  predicted 
temperature  distribution  with  time,  one  inch  above  and  below  the  insiilationj 
for  Sections  "A"  and  "B",  respectively.  Figure  10  indicates  that  the  tem- 
perature below  the  insulation  in  Section  "A"  may  drop  several  degrees  below 
freezing  for  a  period  of  time.  Figure  11  indicates  a  similar  (but  lesser) 
effect  for  Section  "B".  However,  it  is  felt  that  the  designs  are  thermally 
adequate  for  the  following  reasons?  (l)  the  assumed  upper  boundary  conditions 
are  on  the  cor^ervative  side,  (2)  the  assumed  air-surface  transfer  factor 
is  conservative,  (3)  the  lower  boundaary  condition  may  be  on  the  conservative 
side,  and  (k)   temperatures  down  to  28°F  can  be  tolerated  below  the  insulation 
for  several  days  without  damage  (36) « 

A  similar  analysis  was  performed  for  the  Control  and  the  existing  high- 
way section,  the  results  being  shown  on  Figures  12  and  13,  respectively. 
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Figure  12  indicates  that  ftrost  penetration  of  almost  two  feet  into  the 
subgrade  of  the  Control  Sectioa  is  possible.  The  atsalysis  for  the 
Bxistiog  highway  section  was  performed  in  order  to  determine  if  deteriori- 
atioR  of  the  existing  roadway  could  be  partly  attributed  to  frost  action. 
Figure  13  indicates  that  frost  could  Jaave  penetrated  to  depths  of  almost 
three  feet  into  the  subgrade  of  the  SKisting  highway  section,  thus  indica- 
ting that  part  of  the  deterioriation  of  the  present  roadway  may  be  attri- 
buted to  frost  action. 

Thicknesses  of  the  existing  pavement  coo^oneats  vary  from  approximately 
^  to  5  inches  of  bituminous  material  and  approximately  5  to  9  inches  of 
old  road  gravel.  The  thickness  values  \2sed  in  the  analysis  were  median 
values  of  the  abovs,  and  therefore,  the  po£sibi3J.ty  ^cists  that  even  deeper 
frost  penetration  occurred. 

Width  of  Insulation 

In  Section  "A"  the  insulation  will  be  extended  five  feet  beyond  the 
edge  of  pavement,  providing  a  total  insulated  width  of  3^  feet,  while  the 
iasiiLlation  will  be  extended  the  full  width  of  the  shoulder  in  Section  "B", 
providing  a  total  insulated  width  of  U6  feet. 

The  purpose  of  extending  the  insulation  beyond  the  edge  of  pav^nent 
is  to  BKJximiKe  the  volume  of  subgrade  subjected  to  one-dimensional  heat 
flew.  With  reference  to  Figure  k^   this  means  extending  the  insulation  far 
enough  beyond  the  edge  of  pavement  to  greatly  reduce  lateral  heat  transfer 
fr«a  the  subgrade  (the  so-<;alled  "two-dJjEensioHal  effect") .  The  optijnal 
distance  to  extend  the  insulation  to  provide  this  protection  is  beyoM 
aE^lytical  detersvLaation  at  this  time,  Eowever,  current  reeeartih  at 


Purdue  University  is  focused  oa  the  tMo-dimensioaal  problem.  Field  data 
ars  needed  to  verify  and  refine  the  analysis »  With  iiastrumentation  pro- 
posed in  this  report,  interpretations  of  the  heat  flow  patterns  at  the 
sdges  of  the  two  different  irisulated  widths  can  be  made , which  shovild  prove 
(ra3.uat)le  in  the  solution  of  the  two-dimerisional  problem. 
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RECOMMEaiDED  PEBFCOmnCE  EVALUATION  METHODS 
Both  the  thermal  and  structural  performance  of  the  test  sections 
vjill  be  investigated.  If  frost  penetration  through  the  insulation  is 
substantially  prevented,  the  structural  performance  of  the  insulated 
sections  should  prove  superior  to  that  of  the  non-insulated  one.  To 
evaluate  the  thermal  and  structural  performance,  the  follov7ing  methods 
are  recommended. 

Thermal  Performance 

To  measure  temperatures  at  various  strategic  locations  throughout 
the  cross-section  of  each  of  the  three  sections  comprising  the  test  in- 
stallation, thermistors  are  recommended.  Although  thermocouples  have 
been  used  in  similar  experimentation  in  other  locations,  thermistors  are 
favored  for  reasons  of  durability,  stability,  ease  of  taking  measurements, 
and  economy.  The  thermistor  units  will  be  coated  with  epoxy  and  the  leads 
i7ill  be  encased  in  suitable  plastic  tubing  for  protection  against  the  ele- 
ments of  the  in-service  environment. 

The  recommended  number  and  location  of  the  thermistors  for  Sections 
"B",  "C"  and  "A",  are  shown  in  Figures  li+,  15  and  l6,  respectively.  Section 
"B"  contains  39  thermistors;  Section  "C",  2k;   and  Section  "A",  I42...  total- 
ing 105  for  the  entire  project.  As  shoirm  in  the  Figures,  the  thermistors 
V7ill  be  installed  in  vertical  strings  at  the  center  section  of  each  of  the 
three  sections.  It  is  assumed  that  the  thermal  conditions  will  be  symmetrical, 
and  only  one-half  of  any  cross-section  is  instrumented.  As  e:q)lained  in 
greater  detail  later  in  the  report,  all  thermistors  located  in  the  subgrade 
will  be  installed  in  the  side  of  a  fotir  foot  trench,  with  the  exception  of 
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three  thermistors  in  Section  "B",  two  in  Section  "C",  and  three  in  Section 

j 

"A",  These  eight  elements  are  located  at  depths  below  the  bottom  of  the 
trench, and  wilD.  be  installed  in  small  diameter  auger  holes  made  through  the 
oottcan  of  the  trench.  At  each  section,  all  thermistor  strings  will  be  led 
out  beneath  the  shoulder  to  a  terminal  board  on  the  back  slope. 

The  logic  of  the  thermistor  locations  may  be  explained  as  follows o 
Thermistor  columns  A  (Section  "B"),  H  and  I  (Section  "C"),  and  M  and  N 
(Section  "'A")  were  located  to  obtain  icformation  on  the  effect  various 
surface  materials  have  on  the  air-surface  transfer  factor «  Thermistor 
columns  B,  C,  D  and  E  (Section  "B")  and  columns  0,  P,  Q,  R  and  S  (Section  "A") 
were  placed  to  measure  the  two-dimensional  heat  flow  pattern  at  the  edges 
of  the  insulation.  As  stated  previously  in  the  "Design  Procedure"  section, 
it  is  felt  that  valuable  field  data  will  be  generated  by  these  thermistors 
to  help  in  the  development  of  a  two-dimensional  heat  flow  prediction  model., 
Such  a  model  provides  valuable  guidance  in  the  selection  of  width  of  the 
insulating  layer, 

Thermi-stor  columns  G  (Section  "B"),  L  (Section  "C")  and  U( Section  "A") 
were  placed  to  measure  the  one-dijnensional  heat  flow  conditions  at  the  center 
of  the  cross-sections  and  thus  provide  a  check  on  the  predictive  capability 
of  the  extant  one-dimensional  model  (g).  The  readings  ftrom  these  thermistors, 
with  time, can  be  compared  to  prediction  curves  such  as  those  of  Figures  10,  H 
and  12.  The  thermistors  located  immediately  above  and  belw-;  the  insuls-tion 
in  these  columns  are  critical,  and  a  pair  is  provided  to  allow  for  malfuncti'^ns » 
The  two  deepest  thermistors  in  each  of  these  center  line  columns  are  placed 
to  obtain  xmch  needed  information  on  the  lower  boundary  conditions.  DF&a 
from  such  depths  are  largely  unavailable  =  The  importance  of  proper -V  assess- 
ing the  loiter  boundary  conditions  in  the  design  of  an  insulated  pavement  has 
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previously  been  stressed,  with  reference  to  Figure  ?•  The  exact  depths  of 
the  center-line  thermistors,  were  selected  with  the  idea  of  direct  con5>arison 
of  thermistor  readings  between  sections,  as  shown  in  Figure  17 . 

In  Section  "C"  thermistor  columns  J  and  K  provide  information  on  the  heat 
flow  conditions  in  the  shoulder  and  at  the  edge  of  pavement,  respectively. 
In  addition,  when  these  tx*o  columns  are  coupled  with  column  L,  a  general 
picture  of  the  thermal  pattern  for  the  cross-section  is  obtained. 

Thermistor  column  T  (Section  "A")  and  F  (Section  "B")  provide  links 
between  the  temperatvire  pattern  at  the  edge  of  the  insulation  and  the  center 
line  of  their  respective  sections,  and  will  alloi-7  construction  of  thermal 
contours.  In  addition,  column  F  provides  information  on  the  heat  floii?  at 
the  edge  of  pavement. 

Overall,  it  is  felt  that  this  instrumentation  pattern  can  yield  a  su- 
perior description  of  the  thermal  regime  of  the  three  sections,  and  will 
provide  highly  useful  input  for  checking  both  one  and  two-dimensional 
therjmal  prediction  models. 

Structural  PerforEmnee 

It  is  believed  that  for  the  purposes  of  this  project,  Benkeljoan  beam 
tests  offer  the  best  means  of  comparing  the  structural  performance  of  the 
insulated  pavement  sections  with  the  control  section.  Other  means  of  evalu- 
ating structural  performance,  such  as  the  use  of  a  roughometer  or  profilemeter, 
are  concerned  with  determining  a  serviceability  rating,  which  over  the  covixse 
of  15  to  20  years  may  only  vary  from  four  to  two. 

As  reported  by  Turnbull,  Ahlvin,  and  Brown  (31) >  Bsnkelman  deflection 
waasuren^nts  graphically  illustrate  the  occurrence  of  pavement  failures  during 
the  period  of  thaw  in  areas  of  seasonal  frost.  Thus,  Benkelman  deflection 
measurements,  besides  giving  an  accurate  determination  of  the  deflection 
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patterns  of  the  test  installation  sections,  v;ill  allm-;  the  effect  of  seasonal 
variations  to  be  directly  observed.  It  is  recommended  that  these  measure- 
ments be  taken  throughout  the  year... at  the  end  of  construction,  during  the 
fall,  at  the  beginning  of  the  winter,  and,  especially,  several  times  during 
the  spring  warm-up.  The  thermal  instrumentation  will  be  particularly  help- 
ful in  determining  vjben  Benkelman  deflection  measurements  should  be  made. 
Comparative  plots  of  these  measurements  with  tim.e  will  correlate  structural 
performance  with  temperature  and  season. 

Other  Measurements 

Several  other  measurements  are  recommended,  with  one  being  particularly 
important , 

It  is  strongly  recommended  that  a  weather  station  be  set  up  at  about 
the  center  of  the  test  installation  (STA;  103  +  OO)  consisting  of  a  standard 
rain  gage  and  a  7-day  thermograph.  The  thermograph  will  produce  a  continuros 
recording  of  air  temperatures  at  the  site,  which  can  be  interpreted  in  vrri- 
ouB  ways  to  jirovide  needed  input  for  the  upper  boundary  conditions  of  tne 
thermal  analysis.  The  thermograph  should  be  protected  in  a  standard  J«So 
Weather  Bureau  shOlter.  This  shelter  vjhich  protects  the  thermograph  from 
precipitation,  condensation  and  radiation,  has  louvered  sides  to  permit  air 
to  circulate  freely.  The  shelter  should  be  installed  with  the  bo&tom  about 
four  feet  above  the  ground,  and  with  the  door  facing  north. 

Monitoring  of  ground  water  levels  at  the  site  is  also  recommended. 
This  can  be  accomplished  in  a  simple  perforated  pipe  kind  cf  installation 
located  at  a  convenient  point  on  the  lower  backslope.  This  will  provide 
information  on  one  of  the  factors  usually  required  to  have  frost  action  in 
a  subgrade  soil,  namely  a  free  supply  of  water. 
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BSC(MffiNDHD  CONSTRUCTION  PROCEDURES 
General 
Constiuetion  procedures  acd  specifications  should  follow  those  required 
r  the  Indiana  State  Highway  CoBmassion  for  flexible  pavement  roads,  with 
cceptions  to  be  noted  in  this  report  or  deemed  necessary  by  the  Resident 
agineer  at  the  time  of  construction » 

Due  to  the  experimental  nature  of  tM.s  project  and  the  instrumentation 
0  be  installed,  it  will  be  necessary  to  place  special  restrictions  on  the 
iz&  and  type  of  vehicles  and  equipment  alltwed  on  the  test  sections  and 
oatrol  section  during  construction.  Sections  "A"  and  "B"  should  be  closed 
o  all  vehicles  and  equipment  other  than  those  required  for  the  construction 
f  the  sections  for  a  period  from  the  final  preparation  of  the  subgrade 
hrough  the  completion  of  final  paving »  With  the  exception  of  the  area  ad- 
acent  to  the  iristrxunent  locations.  Section  "G"  can  be  open  to  such  traffic 
is  the  Resident  Engineer  would  norioally  permit «  The  area  adjacent  to  the 
nstruiaent  locations  in  Section  "C"  should  be  closed  to  all  traffic  except 
is  required  for  coustruction«  It  will  probably  be  necessary  for  the  con- 
bract  or  to  provide  a  temporary  detour  around  the  test  installation. 

After  paving,  the  staxidard  legal  load  restrictions  will  apply  to  each 
3f  the  three  sections.  In  addition,  the  use  of  lightweight  equipnent  and 
special  load  restrictions  may  be  imposed  by  the  Project  Supervisor  or  Resideat 
Sngineer  for  use  in  the  construction  of  the  test  sections  and  instrumented 
area  of  the  control  section. 

Construction  Procedvires 
A.  Test  Sections  With  Insulating  Layer 
lo   Subgrade  Preparation 

The  compaction  requirements  for  the  soil  on  which  the  insulation 
is  to  be  placed  should  be  those  normally  specified  by  the  engineer  for 
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that  soil.  The  surface  on  which  the  insulation  board  is  to  be  placed 
may  be  leveled  and  smoothed  by  mechanical  equipment  as  directed  by  the 
Engineer  so  that  deviations  from  a  10  foot  straight  edge  are  not  in 
excess  of  plus  or  minus  ^   inch. 

Thermistor  Installation  (Subgrade) 

Hand  labor  and  equipment  for  excavating,  backfilling  and  recom- 
pacting  a  four  foot  deep  trench  (needed  for  installing  the  instrumenta- 
tion) across  one  half  of  the  roadway  of  each  of  the  test  sec<,ions, 
should  be  supplied  by  the  contractor.  The  trench  to  be  excavated  in 
the  subgrade  should  be  backfilled  with  excavated  subgrade  material, 
and  recompacted  to  meet  the  normal  specification. 

Installation  of  the  instrumentation  is  not  the  responsibility  of  the 
contractor,  but  will  be  performed  by  staff  ftoa  the  Research  and  Training 
Center  of  the  Indiana  State  Highway  Commission,  West  Lafayette,  or  the 
Joint  Highxijay  Research  Project.  With  reference  to  Figures  ik,   15,  and  l6, 
it  will  be  noted  that  in  thermistor  columns  G,  L  and  U,  the  ti7o  deepest 
thermistors  are  located  at  depths  below  the  bottom  of  the  excavated  trench, 
while  columns  F  and  T  also  have  their  deepest  thermistor  belov?  the  bottom 
of  the  trench.  These  thermistors  will  be  installed  by  loii/ering  them  to 
their  specified  locations  through  small  diameter  holes  angered  in  the 
bottom  of  the  trench.  The  subgrade  material  removed  from  the  holes  will 
be  replaced  and  recompacted.  These  holes  will  be  mavie  by  state  forces, 
using  state  equipment  (either  hand  or  power  auger). 

A  period  of  several  days  should  be  stiff icient  to  excavate  the  trench, 
auger  the  holes,  install  the  thermistors,  backfill  and  recompact.  Plastic 
coverings  should  be  available  to  protect  the  trench,  if  rain  occurs  or 
threatens . 
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3o   Placing  Insulating  Board 

After  the  thermistors  are  placed  ia  the  subgrade,  the  contractor 
will  provide  the  labor  and  eqid.pment  to  install  the  insiilating  material 
upon  the  subgrade  as  directed  by  the  Beaident  Engineer »  The  insulation 
will  probably  be  supplied  in  sheets  approxinstely  2'  x  8'  x  1"  (Section 
"A")  and  2'  X  8»  X  34"  (Section  "B"),  The  boards  will  be  butted  together 
and  fixed  in  position  with  a  minimum  of  tt^o  wooden  skewers  (approximately 
6"  by  3/8"*)  psr  sheet o  The  skewers  should  be  driven  through  the  insu- 
lation into  the  soil  at  an  angle  to  force  the  board  against  the  joint 
being  formed, and  until  flush  with  the  surface  of  the  previously  placed 
boards o  Plgure  l8  illustrates  thiSo 

Placement  of  the  insulatiiog  boards  should  begin  at  one  end  of  the 
test  section  and  proceed  down  the  length  of  the  section «  The  first 
row  of  boards  is  placed  down  the  centerline  of  the  roadway,  using  a 
striEgline  to  insvire  straight  alignment  of  the  boards.  Placement  of 
the  insulating  boards  should  proceed  from  the  center  line  outward, 
with  the  first  row  of  boards  always  remaining  ahead  of  the  outer  rows 
as  placement  advances.  This  is  illustrated  on  Figure  I9,  In  addition, 
all  transverse  joints  should  be  staggered »  ae  shown  in  Figure  20  <.  To 
accoiaplish  this,  each  adjacent  row  will  be  alternately  displaced  one- 
half  the  length  of  the  boards  (k   feet)  for  the  full  insulated  width. 
Figure  21  illustrates  this  fingered  pattern. 

The  area  to  be  covered  is  200"  by  ±  S^f'  (68OO  square  feet)  for 
Section  "A",  and  200"  by  ±  1*6"  (9200  square  feet)  for  Section  "B". 
Previous  experience  indicates  that  the  insulated  board  can  be  placed 
at  about  the  rate  of  800  bd«  ft.  per  man  per  hour  (36). 
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k »   Placenent  of  tfaterial  on  the  Insulating  Layer 

As  soon  as  the  insulation  is  installed,  the  contractor  should 

start  construction  of  the  overlying  base  or  subbase.  These  courses 

will  be  constructed  in  the  following  manners 
Section  "A" 


The  cc»nj>acted  thickness  of  Type  "I  subbase  will  be  not  less  than  6 
inches.  Accordingly,  about  an  8-inch  loose  lift  should  be  placed  and 
compacted.  Coarse  aggregate  does  not  cause  any  appreciable  damage  to  the 
surface  of  the  insulatiOKj  as  long  as  the  top  size  of  the  material  is  less 
than  approximately  four  inches  and  the  material  is  well  graded  (36). 
Additional  pertinent  details  of  construction  for  this  first  lifb  are  as 
follows.   The  material  should  first  be  end  dumped  adjacent  to  but  not 
ugoa  the  insulation  board,  as  shown  in  Figures  22  and  23 <,  It  should  then 
be  pushed  onto  the  boards  and  spread  by  a  lightweight  track  vehicle  such 
as  a  front  end  loader.  The  vehicle  should  be  equipped  with  street  pads 
and  operated  in  a  manner  such  that  at  no  time  will  it  rest  directly  upon 
the  insulation,  see  Figures  2^  and  25.  Placement  of  material  will  be  f-.'cmi 
one  end  only  arid  proceed  in  the  same  direction  the  entire  length  of  the 
section.  Placing  material  from  both  ends  simultaneously,  may  cause  the 
boards  in  the  center  of  the  section  to  buckle  up.  Cosipaction  of  granular 
material  on  the  insulation  should  be  by  equipment  which  exerts  no  more 
than  an  80  psi  contact  pressure  (37).  Coni>action  of  granular  material  on 
the  insulation  should  be  extended  to  the  full  width  of  the  iasulation 
before  subsequent  operations  are  performed,  and  should  meet  the  normal 
specification  requirements.  There  are  d£  special  restrictions  on  the 

3<»  Extreme  caution  should  be  exercised  by  all  equipment  operators  to 
avoid  damage  to  the  insulation  and  instruments. 
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placement  of  the  9-iach  base  which  overlies  the  subbase.  Conventional 
equipment  required  to  complete  the  constaTuction  of  the  section  may  be 
allowed  to  operate  on  the  base,  once  compacted o 

Section  "B" 


A  6  inch  thick  compacted  lift  (or  8  inches  thick  loose)  of  Type  "P" 
compacted  aggregate  base  -rfill  be  placed  and  compacted  upon  the  insulation. 
A  3  inch  thick  (compacted)  lift  of  the  same  material  should  then  be  placed 
and  compacted.  The  construction  methods  for  the  first  6  inch  compacted 
lift  will  be  the  same  as  those  required  for  Section  "A"  in  the  preceding 
paragraph.  The  3  inch  compacted  lift  may  be  placed  in  a  normal  manner. 
Alternately  the  entire  9  inch  base  may  be  placed  and  compacted  as  a  single 
lift,  except  that  the  special  procedures  described  for  Section  "A"  will 
apply... since  this  is  the  first  lift  above  the  insulation. 

5.  Thermistor  Installation  (Base  Course  To  Wearing  Course) 

A  trench  will  not  be  required  to  install  thermistors  above  the  insu- 
lating layer.  As  each  such  layer  is  completed,  construction  of  the  next 
layer  will  be  deferred  until  state  personnel  have  placed  the  instrumentation 
in  the  newly  completed  one.  Placement  of  the  instrumentation  can  be  through 
small  holes  (large  enough  for  a  man's  hand  and  lo\jer  arm)  made  at  the  ther- 
mistor locations  in  each  layer.  These  holes  will  be  refilled  with  material 
similar  to  that  taken  out. 

6.  Bituminous  Base  Course  and  Surface  Course 

The  Resident  Engineer  will  need  to  malce  certain  value  judgements  in 
the  use  of  heavy  equipment  to  compact  these  layers.  The  principal  item 
of  concern  will  be  damage  to  the  thermistors  and  leads. 
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B.  Control  Section 
The  construction  schedule  for  this  section  will  be  similar  to  the 
isulated  ones,  except  that  the  standard  pavement  section  will  be  used. 
le  time  relegated  to  the  instrumentation  of  the  various  layers  will 
tout  the  same  as  that  for  the  insulated  sections.  Normal  construction 
sthods  and  equipment  will  be  used,  except  that  the  Engineer  should  be 
.ert  to  activities  which  could  damage  the  thermal  instxtuaentation. 
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FURTHER  RECOMMENDATIONS 

lo      At  the  time  of  construction,  or  before,  representative  samples 
of  generous  proportions  should  be  taken  of s  (a)  each  soil  type  used  in 
the  sections  as  conrpacted  subgrade;  (b)  each  soil  type  loaded  in  situ 
as  pavement  foiindation,  to  the  depth  of  the  deepest  thermistor;  and 
(c)  each  granular  material  used  as  base  or  subbase. 

2o      Selected  samples  obtained  in  (l)  above  should  be  subjected  to 
a  program  of  laboratory  evaluation  for  determining;  (a)  standard  indices 
of  classification,  compaction,  and  load-deformation;  and  (b)  such  thermal 
characteristics  as  are  deemed  practicable  and  desirable. 

3.   A  program  of  data  collection,  reduction  and  analysis  should  be 
established  which  will  serves  (a)  to  correlate  the  pavement  sections 
performance  with  their  traasieEit  environment;  and  (b)  to  validate,  refine, 
and  extend  analytical  models  for  thermal  pavement  design, 

U.   A  detailed  program  for  implementing  the  recommendatioas  of  (1,2,3) 
above  should  be  developed  jointly  by  the  Research  and  Training  Center  and 
the  Joint  Highway  Research  Project,  and  be  submitted  to  the  appropriate 
State  authority  for  review  and  action.  Since  construction  of  the  test  facil- 
ity is  imminent,  there  is  soms  lix-gency  to  this  recommendation.. 
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I-««««»»»»»*»»*#»»»»»«#»«»***#»«»*«**«#«»*«»»»«.«»»»»«#»1HHHHJ«»»»»»«#«» 

^♦      THIS  PROGRAM  is  PREPARED  FOR  THE  PREDICTION  OF  TEMPERATURE   «» 
DISTRIBUTION  IN  AN  ARBITRARY  MEDIA  UNDER  VARI0U5.  INITIAL  AND    »» 


h« 


**      BOUNDARY   CONDITIONS      -=ONE    DIMENSIONAL   ANALYSIS—  «♦ 

►«»»«»♦*«»»«♦♦♦♦»♦♦♦*»*«♦*»♦*#»»♦«#»#»»«»»»♦«»*»«»»««♦««««-»»♦#»»»*«»» 

DIMENSION    Ta6,500)»IOC90J 

DIMENSION    ZC15J  ^OENUS  J  »WC«  15  )  sZZUS  )  »SPHT«  15  )  »C0NDa5  > 
DIMENSION    Cia5J»C2(15j9C3C15S»C4(15)»C5«153,C6il5)»C7a5) 
DIMENSION    STEPAa5)»STEP8C15USTEPTa5)  eTFUS  J  e  KE  ( 15  J 
DIMENSION    STEBAS999}  .STEBB(1999  )  ♦STEBT  (999  )  9TB(999  }  ,TBLC999  J 
COMMON   Cl»C2»C3f.C4»C5sC6.C7»T,SpHT»C0NDi.DENi.WC,J»K,M»Z»D 

1  »TB»TBL»DTFOAM»JF0AM,LB0ND 

DATA   CGS»FPH,CENT  »FAHR/3HCGS  93HFPH,4HCENT  i»4HFAHR/  »TFREZ/0./ 

1       .CELL/SOO./ 


**     SOIL  PROFILE  AND  ITS  PROPERTIES   *«*««♦**«**««*««#*««»«»»♦♦»»♦♦«*» 
►♦    ONE  DIMENSIONAL  ANALYSIS        ♦«»*♦#»»#»«♦»#«#*«♦*«»*«»«««♦»♦♦♦ 


SUBSCRIPTS  Je  K 
JtK^X  »2»3«»** 


DENOTING  SPACE  AND  TIME  STEPS  RESPECTIVELY 


)0  READ(5t9>  ID 

READ«5,10)  UN!TS»TICE,TEMPI »TEMPB, TEMPO 
READC5»11}  M 

READ«5,12J  «ZUUDEN«J}»WC«JJ»  J«1»M? 

CHECK  STOP  FOR  NEGATIVE  OR  ZERO  THICKNESS  AND  DRY  UNIT  WEIGHT  AND 
NEGATIVE  WATER  CONTENT 

DO  70  J«1»M 

IF<Z<JJ.LE.0.0)  GO  TO  71 

IFCDEN(J).LE.O.O)  GO  TO  72 

IF(WC(J>.LT.O.OJ  GO  TO  73 
to    CONTINUE 

GO  TO  7h 
n   WRITE(6»755 
15    F0RMATUH1,39HNEGATIVE  OR  ZERO  THICKNESS^  CHECK  INPUTI 

STOP 
12    WRITE(6»76J 
76  F0RMATaHl,45HNEGATIVE  OR  ZERO  DRY  UNIT  WEIGHT*  CHECK  INPUT) 

STOP 


60 

73  WRITE<6»77) 

77  FORMATaHl,35HNEGATlVE  WATER  CONTENT.  CHECK  INf*UT) 
STOP 

74  WRITE(6.101J 

103  FORMAT! 1H1965HDATA  FOR  PREDICTING  OF  FROST  PENETRATION  INTO  A  SOIL 
l»WATER  SYSTEM///) 

WRITE(6»9)  ID 

WRITE(6»102)  M 
:02  FORMAT  J //2X30HTOTAL  NUMBER  OF  LAYERS  DIVIDED.  I5///.7X5HLAYER.3X 

1  9HTHICKNESS.3X11HDRY  UNIT  WT.5X13HWATER  CONTENT»10X 

2  20HDEPTH(BEL0W  SURFACE)) 
IF(UNlTS.EQoCGS)  WRITEJ6.103) 
IF<UNITS«EQ.FPH)  WRITEJ6.109) 

L03  FORMAT a9XlHM.9X7HK6/C.M.,10X7HPERCENT.22XlHMS 
LU9  F0RMATa8X4HFEET,9X3HPCF.12x7HPERCENT.2lX4HFEET) 

2ZJl)«0.0 
DO  90  J»1«M 
90  IZ  lJ-¥l )  "IZUi-t-Zi  J  i 


WRITE(6»104)  (J»ZtJ).DEN<J).WC(J),ZZ(J),2Z(J+l) 
04  FORMAT{I10.F13«4,5X.F8.1*8X.F8.2»15XF7«4,3H  =  » 


.  Jal.Mj 
F7.4) 


IF{UNlTS«EQoCGSJ  GO  TO  100 
DO  99  J=1»M 
Z(J)»Z«J)*0,3048 
99  DENCJ)»0EN(J)*16. 01837 


READ  VOLUMETRIC  HEAT  AND  THERMAL  CONDUCTIVITY  OF  ALL  LAYERS. 
WHICH  ARE  ASSUMED.  VARIES  LINEARLY  WITH  TEMPERATURE.  I.E. 

VOLUMETRIC  HEAT  «  C1*(TEMP)  +  C2 

CONDUCTIVITY  =  C3«tTEMPJ  4-  C4 
WHERE  CI,  C2»  C3»  AND  C4  ARE  CONSTANTS 

LOO  READ«5.13)  (CIC J) .C2{ J) .C3|J) .C4( J ) »  J«1»M) 

WRITEC6.105) 
.05  FORMAT </////7X5HLAYER.9Xl5HVOLUMETRIC  HEAT.10X12HCONDUCTI VITY/) 

IF(UNITS.EO.CGS)  WRITE{6»115) 

IFCUNITS.EQ.FPH)  WRITE(6»116) 
l15  FORMAT* 23X1 1HKCAL/C.M.-C»13X11HKCAL/M°HR»C) 
IK    FORMAT  5 23X10HBTU/C,F.°F, 13X1 IHBTU/FT-HR-F) 

DO  110  J*1,M 
110  WRITE(6.114>  J9C1(J).C2IJ)  »C3<JMC4(J) 
L14  F0RMATCI10,10X2HC«»F7.2»4H  T  +.F7,3»4X2HK«.F7o2»4H  T  +»F7.3) 

IFCUNlTS.NEeFPH)  GO  TO  107 
DO  98  J-l.M 

C2(J)=CC2(J)+CltJ)«32.0)«16«05 
Cl( J)«Cl(J)«16e05 
C4CJ)«CC4(J)+C3(J)»32,0)*1.488 
98  C3(J)sC3CJS»1.488 

♦♦  ICE  FORMATION  CHARACTERISTICS  »♦» 

♦»    ASSUME  THAT  FOR  EACH  LAYER.  THE  PERCENT  MOISTURE  FROZEN  IS    »*» 


•*  AN  exPCMANTIAL  FUNCTION  OF  TEMPERATURE,                        ••• 

*«  PERCENT  MOISTURE  FROZEN  -  C5  »  EXP(  C6  »  TEMP  -  C7  »        til 

**  WHERE  C5,  C6,  AND  C7  ARE  CONSTANTS                 ^  *        !!! 

*♦  READ  1  FOR  KNOWN  COEFFICIENTS                                  * 

*♦         I    F§R  SRyM?ER^  •'•^^^"^  ^^^°  ^'^  EXPONANTIAL  FUNCTION    »♦» 
.♦♦♦♦.**>»»*»*,»*,»«»,»»»^,,»^,,.,^,,,^,,^^^^^^^^^^^^^^^^^^^^^^^^»«^ 

107  DO  '08  J«1,M 

RFA0(5,11)  ICE{J» 

i:ej»icecj) 

^0  TO  C121»122»123},  KEJ 

I2r  IFtTICE.EO.CENT)  IC0RF«1 
IFCTICE.EQ.FAHR)  ICORF-2 

CALL  CURVE(C55,C66,C77,IC0RF) 

C5|JJ«C55 

C6(JJ«C66 

C7CJ)«C77 

GO  TO  108 

21  READ(5,12)  C5 C J) f C6t J) ,C7( J) 

GO  TO  108 

23  C5CJ1I*1.0 
C6( J)«0,0 
C7(J)aO.O 
GO  TO  108 

08  CONTINUE 

WRITE(6»117) 

^'',l^^^'^Ji^^^^^'7>^^»^^y^^*exU3HlCE    FORMATION  CURVE  IN  EXPONANTIAL  FUNC 
IT  ION  I 

IF(TICE.EQ.CENT)  WRITEC6»125» 

IFCTICE.EO.FAHR)  WR!TE«6.126 ) 

25  FORMATC33X20HTCJ.K)  IN  CENTIGRADE) 

26  F0RMATC33X20HTCJ.K)  IN  FAHRENHEIT) 
DO  106  J«1,M 

IF(ICE<J)«EQ.3)  50  TO  124 
WRITEC6.1181  J,C5<JS,C6«J),C7(J» 

X8  FORMAT?  no, 8X19HPERCENT  WC  FROZEN  «,F6.2,7H  -  EXP(,F6.3.Uh  «T(J»K 
1)  ♦  »F6,3»2H  )) 
GO  TO  106 

24  WRITE(6»119J  J 

19  FORMATCI10,8X62HNO  ICE  WILL  BE  FORMED  IN  THIS  LAYER,  SINCE  NO  WATE 
IR  IS  PRESENT) 
GO  TO  106 
06  CONTINUE 

IF(TICE.EO.CENT)  GO  TO  200 
DO  95  J»ltM 

C7(J)«C7(J)+C6(J)«32.0 
9*    C6(J)«C6CJ)»9. 0/5.0 

»*  INITIAL  CONDITION  fTCJ»l)     J»l»2f.»M  •• 
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*♦    READ  1  FOR  CONSTANT  «» 

*♦         2  FOR  ALGEBRAIC  FUNCTION  ♦# 

«♦        3  FOR  TRIGONOMETRIC  FUNCTION  ♦♦ 

«♦        h    FOR  STEP  FUNCTION  »» 

**         5  FOR  DATA  TO  BE  FITTING  INTO  AN  ALGBRAIC  FUNCTION        ♦» 

*♦         6  FOR  KNOWN  TEMP.  AT  MID^POINT  OF  ALL  LAYERS  •• 

200  READ«5»11)  INIT 

GO  TO  «201,2Q2t203»20<»»205,206).  INIT 


I   T<Z.1)«C0NSTANT  FOR  ALL  2 


IQl    READ(5»12}  TINIT 
TIB»TINIT 
DO  221  J«1>M 
iZl    T<J»1)»TINIT 
TBLa)«»TlNIT 
WRITE(6»23)  TINIT»TEMPI 
23  F0RMATI///2X17H1NITIAL  CONDI  TION»10X»46HI  NITIAL  TEMPERATU^iE  "  CONS 
ITANT  FOR  ALL  DEPTH  «fF10»4»8H  DEGREE  .A^) 
GO  TO  299 


I   T(Z»l)«A"»-B*Z-t-C»2»*2-»-D»Z»»3 


Oi  READ(5»13y  AlNI»BINIfCiNl»DINI 

TIB»AINI 

DDZ»Z(l)/2,0 

DO  225  J»1,M 

TCJ.l»»AlNI+BlNl*DDZ-fCINl»DDZ»»2+DlNl»DDZ«*3 
125  D0Z*DDZ+Z(J)/2.0+2(J+l)/2«0 

TBLa)sAlNI+BlNl*DDZ+CJNl*DDZ»»2*DlNI*DDZ«»3 

WRITE(6*24)  AINI»BINI*CINI »DINI 
24  F0RMAT(///2X17HINITIAL  CONDI  TI0N»5X7HT  (Z»0)b,F10,4,2H  -♦•»F1o<.4»5H  ♦ 
IZ  +»F10.4»8H  ♦Z»»2  +.F10.4»6H  *Z**3 ) 

GO  TO  299 


I   T(Z»1J*A*SIN«Z*PI0L)  ■¥   B*C0SCZ*PI/L),  I 

I   WHERE  L  IS  ONE  HALF  PERIOD  IN  SPACE  SCALE  I 


03  READ(5«12)  AINI vBlNl pHALFZ 

TIB-BINI 

DDZ«Z(l)/2,0 

DO  229  J«1,M 

2PI*DDZ»3.14159/HALFZ  , 

T{J,l)«AINI»SINeZPI)+6INI«C0S([ZPI) 
29  DDZ-DDZ+Z(J)/2.0+2{J+l)/2.0 

TBLa)«AiNi»siN{ZPn-i-B  INI  ♦cost  zpn 

WRITE(6»25)  AINI.HALF2»BINI>.HALFZ 
25  F0RMAT«///2X»17HINITIAL  CONDI TI0N»5X7HT (Z»0 )«»F10.A.12H  SIN(3,i4»Z 
l/F5#2t3HJ  •f»F10.4»12H  C0S<  3.  14«Z/F5.2.1H)  ») 
GO  TO  299 
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STEP  iNlTUL  CONDITION   T{ZfI}»A.  Z«ZO  TO  Z«21,       T(Z»1|-B.    I 
Z-Zl  TO  Z-Z2.  •••   T(Ztl)«N»  Z»Z(N-1)  TO  Z«ZN  I 

>A   REA0«5«in    NSTEP 

READ{5,12)    (STEPAdJ.STEPBjDsSTEPTd),    I»1,NSTEP> 

TlB«STEPrClJ 

DDZ«ZCl)/2,0 

I-l 

DO   232   J»1,M 

»3    IF«DOZ.GE.STEPA«n.AND.DDZ»LT,STEPB(I))    60    TO    23A 

I-I+l 

GO   TO   23: 

14    T(JtlJ»SlEPT(n 

12    ODZ«DDZ+2U)/2oO+2(J+l)/2.0 
1  =  1 

(5    IF(DDZ.GEBSTEPA(I).AND.D02.LT,STEPB(in    GO   TO   236 

Jal+l 

GO   TO   235 
16    TBLdJ-STJPTU) 

WRITEJ6»2".) 
!7    F0RMATJ///?X42HINITIAL   CONDITION    IN    FORM    OF    STEP    FUNCTION) 

WRITEJ6»28)  CSTEPAcnBSTEPBcn,STEPT«I),       I«1,NSTEP) 
18    FORMAT (/7X"HFR0M    Z»fF10.4,6H    TO   Z»,Fl0.4,2xl3HTEMPERATURE    «»F10.4) 

GO    TO   299 


INITIAL    CONDITION    IS   A    SET    OF    DISCRETE    NUMBERS 


•    CALL    FIT(AIH,8lNI,ClNI.DJNn 

TIBaAlNl 

DDZ»Z(1)/2,Q 

DO    245    J«ls1 

T<J»II«AINI9-BINI*DDZ+CINI*DDZ*»2+DINI«DDZ*»3 
>    DDZ«D0Z+Z(J)/2,0+Z(J+l)/2,0 

T8L ( 1 > » A I N I +B I N I ♦DDZ+C I N I *DDZ*»2+D I N I ♦DDZ**3 

WRITE(6«29>    AINI.BINI.CINI ,DINI 

F0RMATC///'X46HINITIAL  CONDITION    BY    CURVE    FITTING    AS      T(Z.O)' 
1      F10»4»2H  +,F10.4,5H   »Z    +bF10«4»8H    «Z**2    +,FlO«4,6H   ♦ZW^SJ 

GO    TO   299 


TEMPERATls^ES  AT  CENTERS  OF  ALL  SUBLAYERS  GIVEN,  (VALUES  MAY  BE    I 
OBTAINED  FROM  PREVIOUS  COMPUTATION)  I 

READ(5,1  )  TIB,(T(J»1),  J«l ,M) .T8L( 1 ) 

WRITE(6fiO)  TIB»  (T(J»1)»  J»1,M) 
»0  F0RMAT(///2X17HINITIAL  CONDI TION/5X102HTEMPERATURES  AT  CENTERS  OF 
lALL  SUBLAYERS(DZ)  GIVEN,  (VALUES  MAY  BE  OBTAINED  FROM  PREVIOUS  COM 
2PUTATIONJ/15F6.1) 

GO  TO  29 » 

IF( TEMPI. -Q. CENT)  GO  TO  298 
TlB«(TIB-32.0)»5, 0/9,0 
DO  297  Jsl,M 


6k 

■I    TIJ»l>«CT«J»l)-32«0J»5«0/9«0 
TBL CI »«<TBL(l)-32.0)»5. 0/9.0 

»♦  TIME  INCREMENT.  DURATION  OF  FREEZING  PERIOD  ♦»♦ 

»8  READ(5,13)  THOUR.DT »DTOUT, TOTAL 
REA0(59l5J  JFOAM,DTFOAM 
WRITE{6»31S  DT 

Jl  FORMATU////2X23HINCREMENT  OF  TIME*  DT  *»  F7.3.7H   HOURS) 
ITIME«DT*CELL/2*.0 
TIME*ITIME*24 

IF«TIME.LTe24.0>  TIME»DT»CELL 
TOUR«THOUR»TOTAL 
IFCTlMEoGE.TOURJ  TIME»TOUR 
LTIME»T0UR/TIME=»-0«99 
NDAY«TOUR/24.0 
TOAY««NDAY 
NH0UR^T0UR-TDAY*24, 0 

WRITEC6»32]|  TOURsNDAYeNHOUR 
J2  F0RMAT<2X24HT0TAL  DURATION  OF  TIME  »»F7«1,8H  HOURS  ••I-'^* 
1   9H  DAYS  AND.  I3»6H  HOURS  J 

KT»TIME/DT-*1.0 
KTKsKT-l 

t*  BOUNDARY  CONDITIONS,  TCOgK)  OR  TB«KJ   K.«l,2,,.,  »♦ 

>«  READ  1  FOR  CONSTANT  *• 

»  2  FOR  ALGEBRAIC  FUNCTION  *« 

^»  3  FOR  TRIGONOMETRIC  FUNCTION  »» 

^»  4  FOR  STEP  FUNCTION  «« 

^*  5  FOR  DATA  TO  BE  FITTING  INTO  AN  ALGBRAIC  FUNCTION        »» 


REA0«5,11»  IBOND 

GO  TO  «30l93029303»304»305),  IBOND 


TC0»TIME»«C0NSTANT  FOR  ALL  THE  TIME 


»1  READ«59l2)  BOND 

DO  321  K*leKT 
'1  TBfKJ^BOND 

WRITE<6»33)  BOND,TEMPB 
'3  FORMAT  «//'lXl9HB0UNDARY  CONDI  TIONSs9X53HTEMPERATURE  AT  BOUNDARY  -  C 
lONSTANT  FOR  ALL  THE  TIME  a^F10.4»8H  DEGREE  »A4) 

60  TO  399 


T  ( 0  »  T I  ME  )  ■  A+B»T  IME4-C*  T I  ME»»2+0«Tl  ME»«3 


I 


6*5 
DDT-START  »o»i^i  i  »o  imr  i 

DO  324  K«1,KT 


TB<K;«ABOND-»-BBOND»ODT/UNlT*CBOND»(DDT/lJNlT)«»2+DBOND*{ODT/UNlT)»*3 
|!^  DDT«DDT+DT 

WRTTE(6f34)  ABOND.BBOND.CBONDsOBOND 
^FORMAT  8 //1X19HB0UNDARY  CONDI TIONS,5X10HT (0,TIME! « .F10.4»2H  +.F10.4 
lf3H  ♦TIME  +»F10,4.11H  ♦TIME»*2  '»tF10.4t9H  ♦TIME»»3' 
GO  TO  399 


T(0,TIME)»A»SIN<PI«KT/T)  +  8»C0.S  JPI«KT/T  J  »  I 

WHERE  T  IS  ONE  HALF  PERiOD  IN  TIME  SCALEf  I 

KTIS  THE  TIME  AT  INTERVAL  K  I 

>  READ«5»14)  ABOND,BBONDfUNIT»STARTfHALFT 
DDT-START 
DO  328  K«i,ICT 
TP I «3. 14 1 59» «  DDT/UN I T  J /HALFT 

TBtK»«ABOND*SIN<TPI)+BBON0*COS(TPI) 
J  DDT-DDT+DT 

WRITEC6»35)  Ab6nD, HALFT »BBOND, HALFT 

>  F0RMAT(///lX»i8HB0NDARY  CONDITIONS. 5X10HT  {0,TIME>-.,F10*4,12H  SINJ3 
1.14»T/F5.2,3H)  ♦,F10,4.12H  COSr3.l4«T/F5.2. IH) » ) 

60  TO  399 


STEP   BOUNDARY   CONDITION       T«0.riivi£)*At    T-TO    TO    T«T1«       T(O.TIME)«BtI 
T«T1    TO    T«T2 T«O.TlMt|»N»    T»T(N-1>    TO    T«TN  I 


y   READCS.lli    NSTEB 

REA0(5»12)     tSTEBACn.STEBBjn.aTEBTd),    I^l.NSTEBJ 

READ«5»12J  UNIT.START 

DO  330  I»l, NSTEB 

STEBA ( n -STEBA ( I ) »UNI T 
)  STEBB(I)»STEBB(I)»UNIT 

WRITE(6»37) 

F0RMAT8///1X,44HB0UNDARY  CONDITIONS  IN  FORM  OF  STEP  FUNCTI0N//5X» 
I  lOHNO  OF  STEP99X20HTIME  INTERVAL(HOURS) »11X11HTEMPERATURE) 

WRIT6l6»38Hl»STEBAtn,STEB8(I),STE8T(I)  »I«1»NSTEB> 
I  FORMAT(I11.9XF10.3,3H  -  ^F10«3»10XF10.3 ) 

DOT-START 

I«l 

DO  331  K»l9<T 

IFtDDT.GE.STEBAlI J.AND.DDT.LE.STEBBd))  60  TO  333 

I«I+1 

60  TO  332 
\    TB(K1»»STEBTU) 
L  DDT«DDT+DT 

GO  TO  399 


BOUNDARY  CONDITIONS  IS  A  SET  OF  DISCRETE  NUMBERS 


S  READC5.12>  UNIT9START 

CALL  FIT<AB0NDbBB0ND»CB0ND9DB0ND) 
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DOT-START 

DO  344  K*1»ICT 

TB«K)«A80ND+BBOND«DDT/UNIT+CBOND»(DDT/UNIT)**2>DBOND»( DDT/UNIT) ♦•S 
^4  DDT-DDT+DT 

WRITE(6»40J  ABOND,BBOND»CBOND»DBOND 
iO   FORMAT<///1X50HBOUNDARY  CONDITION  BY  CURVE  FITTING  AS   T(OfTIME»«. 
1   F10«4»2H  +.F10,4»8H  ♦TIME  ♦»F10.4tllH  ♦TIME*«2  +,F10.4,9H  ♦TIME* 
2»3> 
GO  TO  399 

►♦  CORRECTION  BOUNDARY  VALUES  FOR  AlR<=GROUND  INTERFACE  ♦• 

>9  READ<5,12)  CORR 

DO  340  IC«1,ICT 

IFITEMPB.EQ.CENT)  TBC K)*TB fK) *9e0/5.0+32«0 

IFtTBCK)  .LE.32.0)  TB«K)'»32.0"=»  32oO=TB  (K>  )»CORR 

IFJTB(K)  .GT.32.0J  TB?)O«32o0+8  TB«KJ-32.0)»(2«0«»CORR  » 
>0  TB(K»«aB«K>-32.0»»5,0/9»0 

WRITEC6»39)  CORR 

>9  FORMAT?/  1X64HC0RRECTI0N  FACTOR  FOR  BOUNDARY  VALUES  FOR  AIR  GROUND 
1  INTERFACE  »♦  F6,3) 

*•   SPECIFY  LOWER  BOUNDARY  CONDITION  ♦♦ 

H»    READ  1  FOR  CONSTANT  «» 

**  2  FOR  PERFECT  INSULATION  •« 

**  3  FOR  TEMPERATURE  SAME  AS  BOUNDARY  LAYER  »* 

►*        4  FOR  SPECIFIED  TEMPERATURE  »» 

>#*»♦»«♦»♦#♦»*»«•#*»♦♦♦**#*#««<»*».»»♦♦»»♦♦»♦#  »^^^^#^^^^.^^♦^^^^#^^^^^^^^^^^»♦»♦«♦««»# 

READ(5tll)  LBOND 

GO  TO  C35l9352»353t354)»  LBOND 

>1  WRITEC6»4i) 

n    FORMATC/lXf62HTEMPERATUREAT  LOWER  BOUNDARY  REMAINS  CONSTANT  AT  AL 
IL  THE  TIMEI 
GO  TO  1110 
>2  WRITE(6t42» 
^2  FORMAT </lXl24HASSUME  THE  LOWER  BOUNDARY  AS  A  PERFECT  INSULATION  TO 

1  THE  HEAT  FLOW»  T{M+1  »)<;»  »TCM»<)  »  WHERE  LAYER  IM+1)  IS  A  FICTITIOUS 

2  LAYER) 
GO  TO  1110 

>3  WRITE<6»43) 

'*3    FORMATS /1X»78HTEMPERATURE  AT  LOWER  BOUNDARY  IS  SAME  AS  TEMPERATURE 
i  OF  M  LAYER,  TBL(K»»TIM»K) J 

GO  TO  1110 
54  REA0<5,12>  BONDL 

WR!TEC6»44)  BONDL »TEMPB 
1»4  F0RMAT</1X»46HTEMP  AT  LOWER  BDY  EQUALS  TO  THE  TEMP  SPECIFIED. 
1   1HUF6.2,1XA4,1H)) 

IFJTEMPB.EQ.FAHR)  BONDL"(BONDL'-32,0)*5. 0/9.0 

GO  TO  1110 

»♦  CALCULATING  TEMPERATURES  AT  VARIOUS  POINTS  ♦♦ 
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LIO    DO   600    LLL«1,LTIME 

IFaLL.EQ.l)    GO   TO    1112 

DO    1114    J»ltM 
114    TCJ,l)«TUfKT) 

TBL(l)«T8LnCTlC) 

REM  A I N  *  THOUR'-TOT  AL 
IF(REMAIN.GT,TIMEJ  GO  TO  370 
TIME-REMAIN 
KT»TIME/OT-»-l,0 
KTK-KT-1 
J7f  DDT«TOTAL 

GO  TO  (311,312.313,314»315|.  IBOND 

Hi    DO   371   K-19KT 
J72    TB(KJ«BONO 
GO   TO   398 

J12  DO  372  K»1,KT 

TB { K ) »AB0ND+9B0ND»DDT/UN I T+CBOND*  C  DDT/UN I T ) ♦♦2+DB0N0» ( DDT/UN I T ) *»3 
J72  DDT«DDT+DT 
GO  TO  398 

113  DO  373  K»1,KT 

TP I«3« 14159* C DDT/UN  IT) /HALFT 

TB|KI«AB0ND»SIN(TPI )+BBOND»COS{TPl ) 
J73  DDT-DDT+DT 

GO  TO  398 

J 14  i«i 

DO  374  K«1,KT 

J84    IF<ODT,GE.STEBA(I).AND.DOT,LE,STEBB(I))    GO   TO   385 
I«I+1 

GO   TO   384 
J85    TB{ICJ«STEBT(n 
J7a   ODT-DDT+DT 

GO   TO   398 

DO    375    K«1»KT 

TB( K  >  »ABOND+BBONO»DDT/UN I T+CBOND* ( DDT/UN  I T ) »»2+DB0ND«  <  DOT/UN I T ) «»3 
J75  DDT-DDT+OT 
GO  TO  398 

J98  DO  397  K»1,ICT 

IF(TEMPBeEQ.CENT)  TB(  K)3TBCIC)»9«0/5oO-t-32.0 

IF(TB(K>oLE.32«0)  TB( K)«32,0=( 32»0»TBU) )«CORR 

IF(TB<K).GT.32.0J  TB(  K)»32.0+<  TB(K)'»32.0)*(  2.0-CORR  ) 
J97  TB(K|»(TB(K)"32.0)»5. 0/9,0 

GO  TO  1113 


I   TEMPERATURE  IMMEDIATELY  AFTER  THE  PROCESS  STARTS  IS  TAKEN  TO      I 
I     BE  THE  AVERAGE   OF  INITIAL  AND  BOUNDARY  TEMPERATURES  I 


112  TB(l>»tTIB+TB(l))/2.0 

WRITE  INITIAL  TEMPERATURES 


•cJ 
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13    WRITEt6p51)    (J»    J»1»M) 

I"0 

DTK-TOTAL 

IF«TEMPO»EQ,CENTS    GO    TO   481 

T IF-TB a i*9« 0/5. 0*32,0 

DO   482    J»1»M 
82    TFcJ)«TUtl)»9,0/5. 0+32^0 

WRITEC6»52)     I .DTK»TlFf (TF ( J)  .    J«1,M) 

GO    TO   483 
81    WRITEt6»52)     l9DTK!.TBtlJ»    «TU»lii»    Jb1»M? 


83    K»l 

80    DO    501    J»1,M 

IFCWC(J).EQ,0.03    GO    TO    502 

IFCT«J9K).6E,TFREZ)    PARTF«0»0 

IF(T(J»K}.LT.TFRE2J    PARTF- «Cf)  f  J)-EXP  K6«  J)»T«  J,K>+C7(  JH  )/l00.0 

GO    TO   503 

02  PARTF«0.0 

03  PARTUal«0«PARTF 


I   COMPUTE  VOLUMETRIC  HEAT 


SPHT(J)«VHrr(J,K),Cl(J),C2«J3,DEN(J),wC<J},PARTF,PARTU} 
I   COMPUTE  THERMAL  CONDUCTIVITY 


C0ND«J)''TC(t(J»K),C3(J).C4CJJ9DENCJ),WCUJ,PARTF,PARTU) 

CHECK  STOP  FOR  NEGATIVE  OR  ZERO  VOLUMETRIC  HEAT  AND  THERMAL 

CONDUCTIVITY 
IFCSPHK JJ.LE.OeOJ  GO  TO  81 
IFCCONDC JJ.LE.OaOJ  GO  TO  82 
GO  TO  501 
31  WRITE<6»83}  J 

53  F0RMAT«////»2X,45HNE6ATIVE  OR  ZERO  VOLUMETRIC  HEATt  CHECK  INPUT, 
I       1HC»12»10H  TH  LAYER) ) 
STOP 

82  wRITE{6«84)  J 

84  FORMATC////»2X»50HNEGATIVE  OR  ZERO  THERMAL  CONDUCri VlTYt  CHECK  INP 
lUTf   lHUI2tlOH  TH  LAYER)) 
STOP 
[)1  CONTINUE 

I   COMPUTE  TEMP  AT  LOWER  BOUNDARY  ACCORDING  TO  SPECIFIED  CONDITION   I 


a  -T"~ll  1  ■  MW-'~-  JO  tl 


IFCK.EQ.l)  GO  TO  411 

GO  TO  (451»452»453.454J»  LBOND 

51  T8LCK)aTBLCl» 
GO  TO  411 


GO    TO   411  ^^ 

TBLCKS-UMtKJ 
GO    TO   411 

TBL<KJ«BONDL 
60  TO  411 


. — =.«.^.<j. 


COMPUTE    TEMP    BY    CALLING    CORRESPONDING   SUBROUTINE    ACCORDING    TO  I 

CALCULATED   DISCRIMINANT  I 

J=l 

IFtJ.EQ.l)    GO    TO   401 

IF( J.Ge«2.0R.J,LT»M5    GO    TO    402 

IFCJ.EQ.M)    GO    TO   403 

T32«(TC2fKlf»Ta»K:)J/(2(2J/COND«2)-fZCl)/CONDas  J 
T21«<T<1»K)=TB(K)  >/«Zf  1 ) /CONDU  >  ) 

SO    TO   404 

:2«C  T(  J-*"!  ,iC)=-T  ( J.O  J  /  i  Z  U+l ) /CONDi  J+:u+Z  <  J) /CONDU » > 
TZleiTCJ.Kj-TCJ^-l.OS/CZCJI/CONDJJJ+ZSJ^D/CONDU^m 
GO    TO   404 

IF(LB0ND.EQ.2)    TBL(K>»T<M»K ) 
T32a«TBL(K)'=T(!M»KM/«ZfM>/CONDCMM 

T2l«(T<M,Kj-T(M='l,Kn/«Z(M)/C0ND(M)+Z(M»-l)/C0N0CM»l)) 
30    TO   404 

[)»T32-T21 

IF« J.EQeJFOAM)    60    TO   406 

IF(D«GT«0,0)    GO    TO   407 

IF(0«LT«0,0>    GO   TO   408 

T(J»K-«-l)«T«J,KJ 

SO    TO   410 

CALL    STYRO«TTJ 
TUfK+lJ^TT 
30    TO   410 

rj<8TCj»ic) 

CALL   WARM(DT»TJK,TT> 

r(J»K+l»«TT 

30   TO   410 

TJK»T{J»K» 

CALL  COOLJDT»TJX,TT) 

r(J,K+ll(»TT 

30  TO  410 

IFJJ.GE-M)  GO  TO  420 

J»J4-1 

SO  TO  400 

»•«#»«««»»«♦♦«♦»»«#♦»♦»♦«»»#»«»♦*»♦«»»*♦*«»♦«»♦•»«♦#♦»»*»♦»»«♦•»♦♦ 

WRITE  OUTPUT  •»♦ 

»•♦♦♦««♦»♦#♦♦»♦«♦«»«♦#»♦«♦♦»#♦♦«♦»»♦«»»*»♦♦»»«»♦«♦♦»»♦♦«»♦♦#♦♦•»♦» 
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10  dtk«dtk;+dt 

ifcdtk.lt.dtoutj  go  to  ^23 
0ut1«dtk/dt0ut 
lout-dtk/otout+0 • 0 1 

0UT2»L0UT 

DIFF«0UT1=0UT2 

IF«ABS<DIFF)<.GT. 0.0001)  60  TO  423 

IF(TEMPO<»E0,CFNT)  GO  TO  421 

TBF»TB  8  K-H)  ♦9,  0/5.0+32,0 

DO   422    J«1*M 
>2    TFU)»TCJ»IC+lJ»9.0/5»0+32.0 

WRITE(6«52J  K»DTK»TBF9  CTFiJ),  J«1»M) 

GO    TO   423 
>1   WRITEC6f52)    K»DTK,TB«  K+1 J  »    (T8J,K-H),    Ja.l»M) 
13    IFCK.6E.KTKJ    GO    TO    601 

K»K+1 

GO  TO  480 
il  F0RMAT«1H1,13H  T/Z  ACC  TIME»3x3HB0Y»2X »15I6 J 
»2  F0RMAT«I4»F8.2»3XF6.2»3X15F6«2) 

)1  TOTAL»TOTAL+TIME 

WRITE«6?53}  TOTAL»DT 
>3  FORMATS /3X23HTHE  TIME  AT  LAST  STEP  »fF8.2»  6H  HOURS. 3Xf 
1   16HCTIME  INCREMENT",  F5,2»7H  HOURS)) 

)0  CONTINUE 


FORMAT  FOR  INPUT 


>:=*<£S>*^  CSS         ^sa^t 


9    FORMAT C18A4 J 

0  FORMAT 8 2XA3»6XA4,6XA4,6XA4,6XA4> 

11  FORMATUS) 

12  FORMAT«3F10.4) 
i3    FORMAT«4F10.4) 

4    F0RMATC6F10»4) 
.5    FORMATCI5»F10,4) 

>«*«««■««««*»«««*»««»«««««'»«•«#'»*«»««*»«*«#«««««•««««»«»•*««*•)(*»»«••« 
GO    TO    800 
END 

■TC    SUBl 

SUBROUTINE    CURVEK»AA,BB»  KORF) 

DIMENSION    XCl5)»Y(15?,Aa00J»8{l00J 

READ<5,11»    N 

P«0.000001 

READ(5,12)     CXCn.YCIJt    I«1»N) 

READ(5,12J    C 

GO    TO    «1»2),    ICORF 

1  A(1J«0,5 
Bfl)«4.0 
GO    TO    3 

2  AC1S«0.25 
B(l)«-4.0 

3  K«l 
JO   F«0. 

G»0  , 
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FA«0. 

FB»0. 

GA«0. 

GBsO. 

DO  101  I»1,N 

E"EXP(A(KJ»XCn+8{K)} 

F»F"f«Y(I  )-C-*'E)*(X(n»E> 

G»G+(YCU"C+E)»e 

FA«FA+X{  I  )*(X<  I  )*YCI  >»E=X(  I  )»E»C+2,0»X  {  n*E»»2) 

FB«F8+X« I )»iY< I J»E=C»E+2.0*E*»2) 

GAeFB 

GB«GB+CYCI J»E«C*E+2«0»E»»2J 

D«FA«GB=-FB»GA 

DELTA««G*FB=F#GBJ/D 

EPS!»{F»GA-=G»FA1(/D 

A«K-»-l)»A5K)+DELTA 

BCK+lJ-BCKj-t-EPSl 

ERRA«ABS « A ( K*l } -A ( K ) J 

ERRB«ABS « B  { K-*- 1 J -8 « K >  J 

IF{ERRA«LE,P.AND,ERRB.LE.PJ  60  TO  102 

K-K+1 

GO  TO  100 

AA=ACKK) 

bb»b«kk:) 
format  u  5  j 

F0RMAT<2F10.4) 

RETURN 

END 

TC  SUB2 

SUBROUTINE  FITCA.BtCfD) 

FITTING  CURVE  BY  LEAST  SQUARE  AND  GAUSSES I EDEL  METHOD 

DIMENSION  SC10).XC365)»Y{365J 
READJSUl)  N 

READC5vl2}  «X«n.Ya3»  I*1»N) 
DO  1  MaltlO 
1  SIMJ^iO, 

DO    2    I«1»N 

sn)«s«i)+Y(n 

S(2|«SC2)-fXCn»Y(  I  ) 
S(3)«S«3»+X«n»»2«Y«n 

sc4)as<*)+x(n»»3»Y(n 

SJ5>«S<5>+xn> 
SC6>»S«6)+X(n»»2 
SC7j«S«7Uxn)«»3 
S<8)«S«8)+XCr>**4 
S«9)»SC9)+x<l )«*5 
S(10)«S(iO)4-X(I  )*»6 

GAUSS-SIEDEL    METHOD   -   SOLVE    LINEAR   SYSTEM 

FN«N 

A«Sa}/FN 

B«SC^J/S«6) 

C«S(3)/S«8) 

D«SUJ/S«10J 


P»0.0001  ^ 

AAbA 

BB«B 

CC«C 

DD=D 

A»SS(l)°B»S(5j-C»SC6>=D«S(tS  > /FN 

B»(SC2)-A»SC5  5'=C»S(7S'=D»S{8n/S(6) 

C»(S(3)~A«S(6)°B«S(7>^D»SC9n/S(8) 

D«CS<4)«=A«SS7)-B»SC8»='C«Sf9!i  J/SUO) 

IF(ABSCA=AA),6ToP)  GO   TO   3 

IFCABS(B-BB)«GT.P)  GO    TO    3 

IF(ABS(C~CC).GT<.PJ  GO    TO    3 

IFCABSCD»DDUGT«P)  GO    TO    3 

FORMAT! 15  I 
FORMAT(2F10.4) 

RETURN 

END 

SUBS 
SUBROUTINE  COOL 5DT»TJK»TT? 

TEMPERATURE  AT  (K+IJTH  TIME  INTERVAL  WILL  BE  LESS  THAN  CKJTH  INT. 

^«#♦♦#«♦♦*»♦**♦*«»#*#«**##»^^*^^#^^^^»*^^»^s.#^t.»^^.»^^^^^^«^^##^^^f^^^^^^^^l^^^^^,^^^^^^,,^^^^ 

DIMENSION  Z«15J»DENC15J»WCa5}.SPHTtl5?»COND(l5)^Ta6»500) 
DIMENSION  Cl«15)»C2C15),C3a5)»aMl5J,C5a5S,C6«:.5j»C7(15J 
DIMENSION  T8(;999J,TBL?9995 
COMMON  CI,C2»C3»C4»C5»C69C7,T»SPHT»C0ND»DEN9WC,J»K.,M»2»D 

»TB . TBL  »0TFOAM  t JFOAM »LBOND 
DATA  TFRE2»P/0..0»0.001/.QL/80«0/ 


IFCJ,EQ«1) 

60 

TO 

110 

IFCJ.LT.M) 

GO 

TO 

120 

IFCJ.EQoM) 

GO 

TO 

130 

UPPER  BOUNDARY  LAYER 


IFCWCCD.EO.O.O)  GO  TO  112 
IFCTJKaLE.TFREZJ  60  TO  117 
DTFaSPHT ( 1 ) ^Z ( 1 ) ♦ ( TFREZ^T JK ) / ( 2.0«D ) 
IFCDTF^DT)  113,111,112 

TT«TFREZ 
RETURN 

A81a2«0*DT/SPHT  C 1 ) /Z  (1 J 

AB2»ZC2J/C0N0f2J 

AB3s2(i)/cON0fl J 

B«A81/AB3 

A»AB1/8AB2+A83) 

TT*A*T«2»K)-(A+B-l,0j»TJK-»-8*TBC>C) 

RETURN 

DT2aDT 

TTaTJK-lO0«0»P 
GO  TO  114 
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13  DT2SDT-DTF 
TJK-TFREZ 
TT«TFREZ~100.0»P 

1'  AB1«2.0»DT2/SPHT(1)/Z{1) 

AB2«2(2?/COND«2J 

AB38«2n.)/CONDa» 

B»A81/AB3 

A«AB1/«AB2+AB3J 

TA8«A*T«2»K)  =  <A+B«1.0)»TJK+B*TB«tC) 

HL«QL»DEN  J 1 » »WC  i I  J /SPHT « 1 ) / lOo .0 

EXPKaEXP(C6{lj»TJK-^C7{lJ  J/ 100,0 
i  EXPKl»EXPtC6Cl)*TT-».C7a  J  J/100,0 

TTT«TT-( TT»TAB+HL«  C  EXPK1  =  £XPK })/<!«  0+HL»C6 1 1 ) *EXPK1 ) 

IFIABS((TTT«=TT),LE.P)  GO  TO  4O 
>  TT«TTT 

GO  TO  115 


INTERMEDIATE  LAYER 


IF«WC<JJ«EQ90.0J  GO  TO  122 
IFCTJK.LE.TFREZ)  GO  TO  127 
DTF«SPHT  ( J ) »2 « J  I » J  TFRE2-T JfC  >  /  (  2 •  0»D ) 
IF«OTF-DTJ  123.121»122 

21  TT«TFRFZ 
RETURN 

ABla2«0*DT/SPHT«J)/Z<J) 

AB2»Z(J+1)/C0NDU+1) 

AB3»2CJJ/COND«J5 

AB4«Z  « J^l i /COND « J-1 5 

B«AB1/«A83+AB4) 

As«ABl/CAB2+AB3} 

TT«A«T«J+lfK)'-«A+B-l,0)*TJK"^B»T(J»l»K) 

RETURN 

n   DT2«DT 

TT»TJK-IOO,0*P 
GO   TO   124 

DT2»DT=DTF 

TJKaTFREZ 

TT«TFRE2=100.0*P 

ABla2.0*DT2/SPHTCJI/2U> 

AB2«Z«J+1J/C0N0CJ+1 J 

A83aZfJ)/COND«J> 

AB^aZ « J-^l  S  /COND  i  J-1 ) 

B»AB1/<AB3+AB45 

A»AB1/«AB2+AB3) 

TAB«A*T «  J+1 »  K  J  "= « A+8~l .  0 )  »T JK+B»T « J-1  »K ) 

HLaOL«DEN ( J  J ♦WC  C  J ) /SPHT  J  J ) / 100 • 0 

EXPRsEXPeCS?  J)»TJK:+C7JJ)  J/100,0 

EXPK;i«exP«C6C  J)»TT*C7C  jn/100,0 

TTT«TT««TT-TAB+HL»(EXPKl-'EXPKj)/«lo0+HL«C6CJ}*EXP<l! 

IFCABS<TTT»TT).LE.P)    GO    TO    40 


"^-T-TTT 
GO  TO  125 


LOWER  BOUNDARY  LAYER 


I  IF«WC{MJ.EQ»0«0>  GO  TO  122 
IFCTJK.LE.TFREZ)  GO  TO  137 
DTF«SPHTCMJ»ZfMJ»«TFREZ=TJK)/C2,0»DJ 
IF5DTF-DTJ  133fl31»l32 

TTbTFREZ 
RETURN 

:  A81«2»0*DT/SPHTIM)/ZCM3 
A82«Z«M)/C0N0«MS 
AB3aZ(M"lS/COND«M-13 
B«A81/«AB2+AB3) 
IF«LB0N0»E0«2}  GO  TO  138 
A-AB1/AB2 
GO  TO  14^ 
A»AB1/C2.0»AB2> 
TBLCKJsTJK 

T  T«A»T8L  S  K )-  C  A4-B-»l  •  0 )  *T  JIC+B«T  { M-1  »K. ) 
RETURN 

DT2aDT 

TT«TJK-100,0«P 
GO  TO  134 

0T2«DT=0TF 

TJK^TFREZ 
TT*TFREZ"100.0»P 

ABla2«0»DT2/SPHTCM3/Z«MJ 

AB2sZ{M)/C0NDfMJ 

AB3«2  «M-1 )  /COND  W-l  J 

B*ABl/jAB2-fAB3J 

IF(LB0ND-EQ.2>    GO   TO    139 

A«A91/AB2 

GO    TO    140 

A*A81/U»0»AB25 

TBL<KlsT«M,|fJ 

TAB«A*TBL  C  K  11  - « A+B-1 .0  5  »TJK*8«T  <  M~l  t K ) 

HL8QL»DEN(MJ»WC{M)/SPHTfMj/100*0 

EXPK«FXP  C  C6  C  M  5  #T JK+C7  i  M  Jl / 1 00 .0 

EXPK1«EXPCC6CMJ«TT+C7((M)  >/100«0 

TTT*TT=(TT==TAB*HL*CEXPKl-'EXPKn/«1.0-«-HL»C6«MJ*EXPKlJ 

IF(ABS«TTT=TT)«LE,PI  GO  TO  4O 

TT»TTT 

GO  TO  135 

I  RETURN 
END 

:  SU84 
SUBROUTINE  WARM  JDT  .TJIC»TT  ) 
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TEMP  AT  tK+l)TH  TIME  INTERVAL  WILL  BE  HIGHER  THAN  tK,TH  TIME  INT. 

DIMENSION  Zn5)»DEN«15).WC(15j,SPHT{15),CONDa5>.T<)6,500) 
DIMENSION  ClC15),C2(15).C3(15),C4{l5),C5(15),C6a5)»C7(15) 
DIMENSION  TBC999J tTBL(999) 

COMMON  Cl•C2.C3.C<^•C5,C6.C7»T,SPHT.COND»DEN»WC,J,^;,M,2»D 
L      »TB»TBL»DTFOAM»JFOAM»LBOND 
DATA  TFREZ,P/0.0,0,001/.OL/8  0,0/ 

IF(J.EQ.l)  GO  TO  110 
IFCJ.LToM)  GO  TO  120 
IF<J.EQ.M>  GO  TO  130 


-«— — 


UPPER  BOUNDARY  LAYER 


IF«WC<l).EQ.O,OJ  GO  TO  117 
IF{TJK.GE.TFREZ)  GO  TO  117 
HL«OL»DEN ( 1 ) ♦WC  C 1 ) /SPHT  CI  J /lOQ.O 
EXPla(C5(l)=EXP5C6Cl}»TJK+C7Cl) J)/100.0 
DTF«SPhT  J 1  »»Z  £  1  >* <  TFREZ-TJK-^HL♦EXP^  )  /  (2,0«D) 
IF(DTF«DT)  112»111.113 

TT=TFRE2 
RETURN 

DT2-DT 

60  TO  118 

DT2»DT=DTF 

TJK«TFREZ 

AB1«2«0*DT2/SPHT(1)/2(1) 

A82«Z<2»/C0ND(2} 

AB3«Z{1I/C0NDC1) 

B«AB1/AB3 

A»AB1/(AB2+AB3J 

TT«A»T(2fK)-«A+B=1.0)»TJK-t-8»TB<K) 

RETURN 

AB1*2.0*DT/SPHT( 1 )/ZCl ) 

AB2«2(2J/C0ND(2» 

A83«Za)/C0ND(l) 

B«A81/AB3 

A«AB1/{AB2+AB3) 

TAB*A*T  (  2  ,K  J  -  ( A-»-B-l  •  0 }  *T  JK+B»TB  C  K ) 

EXPKsEXP(C6«l)»TJIC+C7(l)  )/ 100.0 

TT«TJK=100,0»P 

EXPK1«EXP«C6(1)»TT+C7a)  J/100,0 
TTT«TT=(TT«-TAB-»-HL»(EXPK1-EXPK:))/C1.0+HL*C6C1)«EXPK.1> 

1f(abs«ttt=tt),le.p)  go  to  40 
tt.ttt 

60  TO  115 


INTERMEDIATE  LAYER 

!F<WCU».EO,0,0)  60  TO  127 
IFITJK.GE.TFREZ)  GO  TO  127 
HL«QL*DEN I J  S 'WC ( J ) /SPHT  C  J ) / 1 00. 0 
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EXPJ«i<C5JJ)-EXP(C6(J)«TJK+C7n)3  )/100.0 
DTF-SPHT  ( J  )*Z  C  J  )*  (  TFREZ^T  JK>+HL»EXP J )  /  { 2  . 0»D > 
IF(DTF=DT»    122fl21»123 

TT«TFRE2 
RETURN 

DT2«DT 

GO  TO  128 

DT2«0T-DTF 

TJK«TFREZ 

AB1«2«0»DT2/SPHT(J)/ZJJ) 

AB2«Z « J+1 ) /COND  C  J+1 J 

A83aZU>/C0NDCJJ 

AB4»Z«J~1)/C0NDCJ==1> 

B»ABl/«AB3-»-AB4) 

A»AB1/«AB2+AB3) 

TT«A»T«J+l»|i;)-(A-frB=l,0)»TJIC+B»TCJ-X9K> 

RETURN 

AB1»2.0»DT/SPHT(J>/2IJ) 

AB2aZ  U-»-i  5 /COND « J+l  > 

AB3«2«J)/CONDCJJ 

AB^aZ  C J-1 ) /COND { J°l S 

B«AB1/«AB3+AB4> 

A«AB1/CAB2+AR3J 

TABaA*Tf  J+l»K)  =  (A+B-KO)»TJK-»-B»T«J-l»K) 

EXPIC»EXP«C6C  J)»TJK+C7<J)  )/100,0 

TTaTJK-100.0»P 

EXPIC1«EXP(C6(  J)«TT*C7«J)  J/lOOcO 

TTT»TT»«TT"TAB-^HL*(EXPKl-°EXPK))/aoO+HL*C6CJj*EXPKl) 

IFCABS«TTT«TT).LE.PJ  GO  TO  40 

TT-TTT 

GO  TO  125 


LOWER  BOUNDARY  LAYER 


IF(WCCM!»EQ«0,0)    GO    TO    137 
IFCTJK.GE.TFREZ)    GO    TO    137 
HL«QL*DEN(M)«WC(M)/SPHTCMJ/100.0 
EXPMsCC5CM>=EXPCC6jMl)»TJ»C+C7(l)U/l00«0 
DTF-SPHT  CM)»Z«M)»CTFRE2-TJK4-HL»EXPM  J/ (2.0*0) 
IF(DTF=DT)    132»131»133 

TT»TFREZ 
RETURN 

DT2aDT 
GO    TO    138 
DT2«DT-DTF 
TJIC»TFREZ 

ABla2«0*0T2/SPHTJM)/Z(M) 

AB2«Z<M)/C0ND(M) 

AB3»Z { M-1 ) /COND  C  M~l ) 

B»AB1/(AB2+AB3J 

1FCLB0N0,EQ«2)  GO  TO  139 

A«A81/AB2 


GO  TO  140  77 

A«=AB1/«2.0»AB2J 
TBL«K)»TJK 

TTsA*TBL « K  }«={  A'frB-1,0  J»T  JK+B»T  { M-1  ,>; ) 
RETURN 

AB1«2.0*DT/SPHTCM)/Z«M) 

AB2s2CM;/COND!MJ 

AB3«Z<M-1»/C0ND«M-1 ) 

B*AB1/   A82+AB3) 

IF  JLBC**JD.EQ.2>    GO    TO    141 

A*A8l/AB2 

GO    TO    142 

A»^A8'./(2,0*AB2» 

TBLCOaTJK 

TAB«A*TBL{<)-CA+B=1.0J*TJS<+B«TCM'-lfK) 

EXPK«EXP«C6«MJ*TJ(C+C7«Mn/100.0 

TTvTJK-100.0»P 

F/<PK1»EXP8C6{M)«TT+C7CMH/100,0 

•aT«TT««  TT=TAB+HL*CEXPIC1=EXPK}  >  /(1.0+HL*C6<M)»EXPK1  J 

IF(!ABS(!TTT=TT)oLE6p)    60   TO    40 

TUTTT 

60    TO    135 

RETURN 

END 

SUBS 
SUBROUTINE    STYRO(TTS 

»»*»*«♦**»#*♦*«»«»*«#*#»**«»««•»#»*«♦♦•»*»«««•»«*«#«»»««»»»#»««•»♦«»«« 

COMPUTE    TEMPERATURE    FOR    EXTRA    THiN   LAYER    WHICH   NEEDS   MORE    STEPS   ♦* 

TO   MEET    STABILITY   REQUIREMENT   THAN    OTHER    LAYERS  *♦ 

DIMENSION  ZC15>»DEN<15J»WC«15JfSPHTU5).CONDa5J.Ta6,500J 
DIMENSION  ClC15J,C2a5S.C3Cl5J  tC4Cl5}»C5a5)  »C6(  15  »  sCT  US) 
DIMENSION    TB«999),TBLC999) 

COMMON   ClfC2tC3»C4,C5,C6.C7»T9SpHT9C0ND»DEN»WC.J»K,M»Z»0 
»TB,T8L,DTF0AM»JF0AM>LB0ND 

IF«D.NE..O,0)    GO    TO    10 

TT*T«J.K) 

RETURN 

NDT«DT/DTFOAM+0.0001 

IF( JFOAM.EQ.IJ 
IFC JFOAM.LT.M) 
IFJ JFOAM.EQ.M) 


UPPER    BOUNDARY    LAYER 

Aei»2.0*DTFOAM/SPHT  ID/lili 

AB2a2(2?/COND(2) 

AB3sZ{i}/C0ND(l) 

B-AB1/AB3 

A»AB1/CAB2+AB3^ 

TT«Ttl,K> 


GO  TO 

110 

GO  TO 

120 

GO  TO 

130 
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DO  111  NF0AM«1»NDT 
1  TT»A*T<2»K)-«A+B-l«0)»TT+B*TBl!K) 
RETURN 


INTERMEDIATE  LAYER 

0  AB1»260*DTF0AM/SPHT(JJ/2UJ 
AB2-Z  <  J-H  ) /COND « J+1  > 
AB3»28J)/COND(J) 

AB4»Z  U-1 J /COND  i J^l ) 

B«ABl/(AB3-«-AB4) 

A«AB1/<AB2+AB3) 

TT«TCJ9K> 

DO  121  NFOAMal»NDT 

1  TT«A*T«J-H»ICJ-JA-t-8^1«0)»TT+B«TIJ°l»Kj 
RETURN  ; 


LOWER  BOUNDARY  LAYER 

AB1«2.0»DTFOAM/SPHTCM)/ZCM» 
^B2»Z<M)/C0NDCM) 

AB3s»Z  CM-1 J  /COND«M~l  > 

B«AB1/?A82+A83) 

IFCLB0NDeEQ«2J  GO  TO  132 

A«A81/AB2 

30  TO  133 

A»'AB1/C2«0»AB2J 

TTsTJMsKj 

TBL«K>«T«M»KJ 

00  131  NFOAM^l.NDT 

T T*A«TBL  « K  3 -  ^  A+B-1 »  0  J *TT'»-B»T « M-1  sK  ) 

RETURN 

END 
ITC  SUB6 

FUNCTION  VH«T.A.B9DEN»WC»PARTF.PARTU) 

»   THIS  SUBPROGRAM  IS  WRITTEN  FOR  COMPUTING  VOLUMETRIC  HEAT  OF     »♦ 

!»    THE  LAYER  AT  SPECIFIED  TEMPERATURE  AND  AMOUNT  OF  WATER  FROZEN  *» 
»«»«*#»«#♦«»*♦«»##»**♦»#****««*#»«»»#»««»*«*»»♦»»«»»-**»•«■»»«♦«♦#»»«♦« 

DATA  SPHTW»SPHTI/ioO»0.49/ 

IF«WC.EQ»O.Oj  GO  TO  1 

VH«A*T+B+DEN»CWC/100«0}*JSPHTW*PARTU+SPHTI»PARTFJ 

RETURN 

1  VH«A*T+B 
RETURN 

END 
rC  SU87 
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FUNCTION    TCrr»A»B»DEN»WCfPARTF»PARTUJ 

THIS   SUBPROGRAM    IS   WRITTEN    FOR   COMPUTING    THERMAL    CONDUCT rvlTY      ♦♦ 
OF    THE    LAYER    AT    SPEC8FIED    TEMPERATURE    AND   AMOUNT    OP    WATER  ♦» 

FROZEN  ^^ 

DATA  TCW»TCI/0.514,l,91/»SPWT/2.65/ 

IF(WC.EQ«0,0}  GO  TO  1 

C-DEN/SPWT 

D»DEN*WC/100,0 

TC«C/ « C+D ) » C  A«T+B ) -<-D/  { C+D )  * C  TCI ♦PARTF+TCW^PARTU  J 
RETURN 

TCaA»T+B 
RETURN 

END 


fVJl 


(T 


'OOCOOc;ooa"■?•o■':^T"^7^'^■^x    tctt 

i;OCOOOOOO'^Cr'^C"3~C^7^C    y-  C^   C    T   c 
OOOOC^(7«C7^0C^aoaoop^^^O■•0   r^^^NJ•''^ro(M 

o  o  tr  a-  ■'f-  i'^  ^-  o  i'^  vT  ro  -si  c  c  .-r  i  j-i  ^t  ;r  -~  "M 

'OOC^o^c^cy^■t^c^c^c^c^a^a»aoaococoaoa^cnx 
^i^st-j   -J•N^v}■-J•vJ•>}-sr^t^}•vf^vr-r  j-sj-»i--j- 

•  oco^cr>OLni^r-40'ocvO,}-'\jOCT-r^  o  ~o  jy  .^  -^ 
o  o  T  (T  CT'  J-  a-  o^<roC'CCxixoot>-r-r--*~-*~-t^r^ 


o  no 


Cm    C    — ' 
■.">    INJ    C 


n.   r^ 


•J-    r>  >-sj  -".   X  r^   -c   CO  -J  f^   r-- 


c_ju^>J'cru'U"i)xcyuur-r-r-r^r-.0~yvj>tj>ONO 
^^vJ■^f^J•4•~r^3■^3■J■^~J•^r^T^^r<■>^^<-vT-^ 

oo^in-HXt/^foocr-i^oxvOLOsr-'^o-j-rNjCN, 


o  X  — '  .r  '"vj  f^  -o  >f  ^  ^ 

Cj  —>   Xi  rr^   •T  ^   '-T'  f^  f^J  O 


•^     ^T    -^1    ^    X     .r    vT    1^    ^      3 


CP  cr  r-  "O  -o  in  u'i  in  LO  ^  m  "n  r^  r"  "■-  m  ^  ,^  fT-)  rr^  pr 
ir.  -r>3■^r^3-^>^■sr-J•-r>3■>T^J•■o•Nr^r^r-r■J■r•J• 
^J■'\JO'^x<fL-^J-\^c— ^-C'^o^'^— •'^'-^-'^.T- 
O'^or-r^rox43JDv0'Nif^x— '-opnc<Ni— 'xo 

OX^-Nrm'ncMr\it\j— i^O^H"    i,  O— 'fxi^vj— 'tj— i 

,  lt.  nt  <r  -r  -J"  >3-  <}-  <r  4-  <^  -J-  ro  ,-<•  .^  <r  sj  -r  ^j-  >r  <J"  -j- 

C5Xir\>j-rO-oo>>irvjr\J— •C'CrcriTO'— 'fNifNiOOCJ 
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(ECJTE       IB  JOB 

JJ08 

JFTC  MAIN 

THIS  PROGRAM  COMPUTE  FREEZING  INDEX  FROM  GIVEN  DAILY  TEMPERATURES  *»♦ 

DIMENSION  DAYC365i,TEMP«365)»IO(36»9DD«3651l,DDC8  365UNDAY8  365j 
DIMENSION  xa00»4i»Y(i00»4J,NDATA<4)»XLIN£«20  5»YLINE(20i 
REAL  LOW 
INTEGER  PERIOD 

$00  READC5e95  ID 
READISsslOJ  H 

REA0«5,11J  {DAYdJ.TEHPCn  ♦  Ial»N} 
WRITEC6si9j 
WRITE«6»95  ID 
WRITE(6»20S 

COMPUTE  DEGREE  DAYS  AND  CUMULATING  DEGREE  DAYS 

DO  99  I^IkN 

DDUS»TEMPcn=32.0 
99    NDAYU  )«DAYC  I  S 

ODCaWDDUJ 

DO    100    I»2»N 
00    DOC«n«DOCU-=lJ+DDU) 

WRITE«6»12)     «NDAY<n»TEMP{n  jDDCU^DDCn  J  B    l^lvHi 

FIND    THE    MAXIMUM    VALUE    IN    CUMULATED   DEG    DAY    VS    TIME    CURVE 

HIGHaO, 

DO    200    I«leN 

IF«DDCUUGE«.HIGH}  GO  TO  201 

GO  TO  200 
!01  HIGHsDDCn  J 

IHIGH«I 
too  CONTINUE 

FIND  THE  MINIMUM  VALUE  IN  CUMULATED  DEG  DAY  VS  TIME  CURVE 

LOW«0» 

DO  300  I«1,N 

lF«D0C«n.LT.LOW)  GO  TO  301 

GO    TO   300 
}01    L0W«D0CCKI 

1L0W«I 
)00  CONTINUE 

CALCULATE  FREEZING  INDEX  AND  DURATION  OF  FREEZING  PERIOD 
FI»HIGH==LOW 
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PcRl0D«IL0W~IHl6H+l 
WRITE<6»13)    HIGH»LOW,FI»PERIOD 

9    FORMATUBAA) 

10  F0RMATU5) 

11  FORMAT|10X,2F10.4) 

12  F0RMATjI8»6XF5,l»8xF6.1*i4xF8.1) 

13  FORMAT(///10X31HMAXIMUM  CUMULATED  DEGREE  DAYS  =»F?.l 

1  /10X31HMINIMUM  CUMULATED  DEGREE  DAYS  «»F7.1 

2  //10X16HFREEZING  INDEX  *»F7.1 

3  /29HDURATION  OF  FREEZING  PERIOD  «»I5»///) 

19  FORMAT UHl) 

20  F0RMAT«//5X3HDAY,3X11HTEMPERATURE»3X11HDE6REE  DAYSt 
1   3X21HCUMULATED  DEGREE  DAYS//) 

60  TO  800 
END 
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